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SCREW INSTABILITY OF A PLASMA COLUMN 


F. C. Hoh and B. Lehnert 
Royal Institute of Technology, Stockholm, Sweden 
(Received July 3, 1961) 


In earlier experiments with the positive column 
an instability has been observed in the presence 
of a longitudinal magnetic field.‘~* A theoretical 
explanation of the phenomenon has been given by 
Kadomtsev and Nedospasov* who discuss the 
growth of a screw-shaped disturbance in the plas- 
ma column. Recent experiments support the 
theory.°,® Even if the physical and mathematical 
assumptions of the theory seem to be well founded, 
the mechanism of the instability has not yet been 
fully explained. In particular, the role which the 
ixB force plays is not very clear. In this Letter 
an attempt will be made to interpret the results 
physically. 

According to the theory‘ the plasma is stable 
for a screw disturbance of the form exp(ikz+imy 
-iwt) with m=1 if 








v b 
KX* + FX?+G.+G > 0.163-—-x-£, (1) 
d +r B.D b. 
Oe i 
where 
_1.28+y ,,_0.8(y +2) 
y(y+1)’ _, 
oA ty). 6 effi 
G = G. =0. 1-* or (2) 


and X=kQ1/B,, y= bj(Q7)?/be, Q is the gyrofre- 
quency, 7 the collision time, bg the mobility, v, 
the axial velocity, and Dz the diffusion coefficient 
of electrons. Further, 5; is the mobility of ions 
and 6,a and f,a are the first zeros of J, andJd, 
with a indicating the tube radius. The collision 
frequency of ions is assumed to be much larger 
than their gyrofrequency. 


Examination of relation (1) and the basic equa- 
tions used in the theory* shows that the right-hand 
member contains the imposed electric field and 
represents the effects which destabilize the system. 
The stabilizing effects are provided by the left- 
hand member. In the first term, KX*, are in- 
cluded dissipation effects due to diffusion and con- 
duction of electrons along the magnetic field B,. 
The term Gg represents corresponding effects in 
the transverse direction for electrons and ions, 
and FX? contains a coupling of longitudinal and 
transverse dissipation effects. Finally, the term 
G, arises from a rotation, Eo, Bz /Bz*, of the 
perturbed electron distribution in the unperturbed 
ambipolar field Eo, (see Fig. 1). Since ions are 
assumed to be unaffected by the magnetic field, 
they will lag behind the rotation of the electron 
distribution. 

The sketch given in Fig. 1 shows a density per- 
turbation of screw type superimposed on the un- 
perturbed, steady-state distribution. With the 
directions given in the figure the axial electric 
field E0z tends to “lift up” the electron screw 
relative to the ion screw. This is equivalent to 
a rotation of the electron screw in the positive 
y direction. Due toa subsequent charge separa- 
tion an azimuthal electric field Ey’ will arise 
which tends to drive the particles outwards with 
a speed E y’ xB, /B,? and to destabilize the plasma. 
The pitch of the screw and its angular velocity Oy 
as given in Fig. 1 are consistent with theory* and 
experiments.®»® The outward drift due to charge 
separation resembles that discussed earlier in 
connection with the stability of a plasma bounda- 
ry.’»® In the latter case charge separation is pro- 
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FIG. 1. Left-handed screw instability for m=1. 
The perturbed density distribution of ions is given 
by the screw-shaped body confined by full lines. 
The corresponding electron distribution is indicated 
by dashed lines. 


duced by drift motions across the magnetic field, 
whereas the separation in the present configura- 
tion arises from motions along the field. 





On the other hand, the ambipolar electric field 
tends to rotate the electron screw in the negative 
y direction relative to the ion screw. The field 
Eoy has the direction shown in Fig. 1 under the 
prevailing experimental conditions, i.e., when 
no field reversal takes place near the critical 
magnetic field. This rotation counteracts the 
charge separation produced by the longitudinal 
electric field. In the special case where these 
two effects balance each other the electron screw 
will become steady, with electrons moving in 
screw-shaped paths inside the ion screw. It is 
then also seen that the balance takes place at a 
certain critical strength of the magnetic field. 

In reality the rotation will not be uniform as as- 
sumed in this simple description but the main 
features of the stability situation will still be the 
same. 

Additional stabilizing effects are provided by 
diffusion and conduction which tend to smooth out 
the perturbed density and potential distributions 
and to suppress their asymmetry. 

According to the present interpretation it is not 
unimaginable that the screw instability will exist 
in a large number of situations where there is a 
sufficiently large electric field component along 
the magnetic field. For the mechanism just de- 
scribed, the destabilizing action of this component 
does not depend upon the ionization degree and the 
instability can also be expected to arise in a fully 
ionized plasma. 





1B. Lehnert, Proceedings of the Second United Nations 
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EXPERIMENTAL EVIDENCE FOR BEAM-PLASMA INTERACTION 
IN A LOW-PRESSURE ARGON DISCHARGE* 


Lawrence H. Putnam, f Harry D. Collins, and Norman L. Oleson 
Department of Physics, U. S. Postgraduate School, Monterey, California 
(Received June 30, 1961) 


Akhiezer and Fainberg’ and Bohm and Gross,’ 
in theoretical investigations, have shown that the 
state of an unmodulated beam of electrons moving 
through a homogeneous plasma is unstable. Den- 
sity and velocity fluctuations, arising in the beam 
or the plasma, may excite longitudinal plasma 
oscillations of growing amplitude in the plasma. 
Recent experimental results obtained in this lab- 
oratory tend to corroborate this theory. 

Data have been taken in an argon discharge at 
6.9 microns pressure. The discharge parameters 
were: current 50 ma, tube voltage 40 volts, cath- 
ode-anode separation 40 mm, cathode diameter 
2cm. A cylindrical Langmuir probe was used in 
the discharge at ten positions spaced 4 mm apart 
along the axial length of the tube. Probe potential 
was applied with respect to the anode by using a 
very slow (0.05 cps) triangle wave. Probe current, 
probe current squared, and first and second deriv- 
atives of probe current with respect to probe volt- 
age were plotted vs applied probe potential auto- 
matically by means of an X-Y recorder. Deriva- 
tives were obtained electrically from simple RC 
differentiating circuits where the time constant 
was about 1/100 of the half-period of the driving 
triangle wave. 


Inserts: 
(a) (b) (c) 
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FIG, 1. 


Electron densities and associated oscillations as a function of cathode-probe distance. 


Second derivatives were then transformed into 
energy distribution functions by use of the expres- 
sion® 





3 ad “, 
NOCV)=T\ > 7) Ww” 
p 


where No(eV) is the number of electrons per unit 
of energy per unit volume, i.e., reciprocal elec- 
tron volts and reciprocal cm’, A is the area of 
the probe, mg is the mass of the electron, e is 
the electronic charge, ip is the probe current, 
and Vp is the probe voltage with respect to the 
anode. 

Examination of the distribution functions for the 
ten probe positions between cathode and anode 
shows that two distinct groups of electrons are 
present out to 12 mm from the cathode, a low- 
energy group with an average energy of about 5 
electron volts, and a higher energy group of about 
22 electron volts average energy. At 16 mm from 
the cathode, the distribution function shows that 
the high-energy group of electrons (beam elec- 
trons) is starting to transfer energy to the low- 
energy group (random plasma electrons). All 
other distribution functions at distances greater 
than 16 mm from the cathode show only one group 
of electrons which is essentially Maxwellian in 
form. This fact indicates energy exchange is es- 
sentially complete. Three distribution functions, 
showing the gross changes, are illustrated in in- 
serts (a), (b), and (c) to Fig. 1(A). 
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(A) Solid circles: 


total electron density (from distribution functions); open circles: plasma density (from distribution functions); cros- 
Ses: intermediate energy beam electrons (from distribution function); solid curve: high-energy beam electrons 
{combined technique of J. F. Waymouth |Notes, Massachusetts Institute of Technology Summer School of Plasma 
Physics, 1958 (unpublished); J. Appl. Phys. 30, 1404 (1959)], and G. Medicus (J. Appl. Phys. 29, 903 (1958)]}. 

(B) Ion oscillation intensity associated with curve -x- of Fig. 1(A). 77 
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Densities were computed by numerically inte- 
grating the energy distribution functions to obtain 
the random plasma electron density, beam elec- 
tron density, and total electron density. In order 
to check the magnitudes of density obtained from 
the distribution functions, total electron density 
was also obtained from the slope of the probe cur- 
rent squared vs probe potential curves. Electron 
density as a function of cathode-probe distance is 
shown in Fig. 1(A). 

Frequencies were measured with a cavity wave- 
meter inductively coupled to the probe at anode 
potential. Detected signals from the wavemeter 
crystal were fed into a dc micromicroammeter. 

A frequency of 520 Mc/sec was measured at a 
point 20 mm from the cathode. This frequency 
corresponded with the plasma density, rather 

than the total density, when the relation 


Ww 


W, = p 
0" a-Si 


cited by Kharchenko et al.,’ was applied to it. 
Here w, is the observed frequency in radians/sec, 
w, is the Langmuir plasma frequency, y, is the 
beam velocity, and S is the average thermal plas- 
ma velocity. 

This frequency could only be detected in the vi- 
cinity of 20 mm from the cathode and could be as- 
sociated with a decreasing plasma density gradient 
in accord with the contention by Allis* and the find- 
ing by Mahaffey et al.° that oscillations are initia- 
ted at or near a peak density gradient and grow in 
intensity down a decreasing density gradient. It 


1,2 x 19'0 (Amperes 2) 
£ 3 
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is also significant to note that at this point (20 mm 
from the cathode) the beam density is starting to 
decrease, showing that breakup of the electron 
beam is occurring. 

Another frequency of 472 Mc/sec was observed 
at 24 mm from the cathode. This observed fre- 
quency was also correlated with the plasma den- 
sity and was associated with the decay of another 
higher energy beam of electrons (of an energy ap- 
proximately equal to the tube voltage drop) which 
was present but does not show on the energy dis- 
tribution functions because it was an order of mag- 
nitude smaller in density (see Fig. 1). 

Ion oscillations of about 100 kc/sec were ob- 
served on the oscilloscope throughout the length 
of the tube. These oscillations were greatly dis- 
turbed in magnitude and wave shape in the inter- 
acting region, 16-20 mm from the cathode. A 
plot of ion oscillation intensity as a function of 
cathode-probe distance is shown in Fig. 1(B). 

At a slightly lower pressure (5.9 microns) with 
discharge currents of 100 and 125 ma at a tube 
voltage of 20 volts with essentially the same cath- 
ode-anode separation (42 mm), high-energy elec- 
trons were detected. This phenomenon manifested 
itself in a disturbance to the probe characteristic 
curve at a point in the tube which was associated 
with either a peak density gradient, or a decreas- 
ing density gradient immediately after a peak. The 
disturbances were noted where visual observation 
of the discharge showed that a meniscus was pres- 
ent or where pronounced beam divergence oc- 
curred. This disturbance to the probe character- 
istic-was of two forms and was always evident 


FIG. 2. Probe characteristic 
curve of ion current squared vs 
probe potential, showing high- 
energy electrons which one as- 
sociated with the electron os- 
cillations. 
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well into the ion-collecting region of the curve. 
The first form was a slight change in slope of the 

ip VS Vp Curve at about -100 volts with respect to 
the anode. The slope change became quite pro- 
nounced on the ip’ vs Vp plots. These slope chang- 
es seemed to be associated with a peak density 
gradient, and hence, in accord with the theory of 
Allis* predicting that oscillations are being initia- 
ted at such locations. The second form of disturb- 
ance (only observed at the highest discharge cur- 
rent, 125 ma) was a step increase in probe cur- 
rent at about -50 volts with respect to the anode, 

as well as the slope change at about -100 volts, 

as described above. An example of these disturb- 
ances is shown in Fig. 2. 

Although in this case no frequency measure- 
ments were made to correlate with these disturb- 
ances, from the work of Emeleus and his group 
showing a correlation between oscillations and 
meniscus formation, it appears probable that the 
disturbances are caused by high-energy electrons 
produced as a result of electron oscillations in the 
interacting meniscus region of the discharge in 
which energy exchange is taking place between the 
beam electrons and the random plasma electrons. 

It is suggested that the unmodulated beams of 
electrons, coming from the cathode, passing 
through and maintaining the plasma, suddenly 


become unstable and give rise to electron oscilla- 
tions of the plasma in the presence of an associa- 
ted decreasing plasma density gradient. In this 
interaction of beam and plasma electrons, elec- 
trons having energies of up to about 100 electron 
volts (well in excess of the tube voltage drop) are 
produced. The step-like disturbance at -50 voits 
is interpreted to be a monoenergetic group of 
electrons associated with a minimum of a density 
gradient and hence peak oscillation intensity. 





*Work supported by the Office of Naval Research. 

tCaptain, United States Army. 
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LIFETIME EFFECTS IN CONDENSED FERMION SYSTEMS* 


A. Bardasis and J. R. Schrieffer 


Department of Physics, University of Illinois, Urbana, Illinois 
(Received July 7, 1961) 


It has been pointed out’ that lifetime effects in 
the superconducting state can be of great impor- 
tance in obtaining the excitation spectrum of the 
system. The concept of a clearly defined elemen- 
tary excitation only has significance when the 
state is long lived; this condition is expressible 
by E, «xE, where the energy of the excitation is 
E=E,+iE,. The original* BCS theory of super- 
conductivity obtains exact eigenstates for a re- 
duced Hamiltonian which contains no damping 
effects. If one attempts to include finite-life- 
time effects, a Green’s function formulation of 
the problem provides a more convenient approach 
than the variational method used by BCS. This 
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FIG, 1. Self-energy diagram used to evaluate the 
integral equation for the energy gap Ap: 





technique has been used by many workers, and we 
follow most closely the formulations of Nambu* 
and Eliashberg.* 

Within these schemes, the evaluation of the self- 
energy diagram of Fig. 1 gives an integral equa- 
tion for the energy gap Ap’ at T =0°K as a function 
of the momentum-energy four -vector p = (p, Po): 

A pY p! a) 
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"»'p is the two-body scattering potential, 7=0*, and Ep is the electron kinetic energy measured 
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from the Fermi surface. The quantity Z p is given 
by 1 -(fp/Po); where bp satisfies an integral equa- 
tion similar to (1) and gives the normal state sin- 
gle-particle self-energy in the limit Ap> 0. The 
imaginary part of f» gives rise to damping effects 
in the energy gap equation. If the potential Vp: p 

is taken to be independent of p,’ and p,, the ana- 
lytic properties of the single-particle Green’s 
function G(p) allow us to reduce the integral equa- 
tion to 


dp! odp,’ 4, /Z@,,? 
A 21m fev f te hee —— , (2) 
p (27) p'pJg (21) (by? -E,,P+in 








dp’ 


where we have set 
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E,,) “GP le, P+, F1-€,,/2,". 0 
It is straightforward to generalize (2) to include 
potentials depending on p, and p,’. Assuming that 
the analytic continuation of G(p) into the lower | 
half-plane for p,>0 has a simple pole correspond- 
ing to an elementary excitation of energy 


-~_ FF am 
E , =E(p’,p . pare «6S (4) 
the integral equation becomes 
= 7 
A(p’,E., ) 
(p > 
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As we have noted previously, the damping ef- 
fects are contained in the quantity Z,=1 - (fp/Po)- 
When Z p is set equal to unity, the integrand be- 
comes real and the equation reduces to that of 
BCS. In our calculation, we investigate the effect 
of damping on the integral equation by consider- 
ing various forms of the function ¢,. The real 
part of tp may be incorporated in a renormaliza- 
tion of the real excitation energies €p, so that we 
need only include ImZ,p in Zp: We consider a gen- 
eral form for Z», in which the damping is propor- 
tional to a power of p,: i.e., 

Zz =1 ia lpol” sgn(po) (6) 
p 0 
With this form for Z p? the energy gap equation 
becomes 


A =- ay 
p i PP on 
p 


where po’=Ep ”, We will consider the three cases 
n=0, 1, 2. It should be noted here that our as- 
sumed form for Im, is not entirely satisfactory 
since this quantity must vanish for |po| < Ap. 
Therefore, in the n=0 case we may only consider 
the Ap> 0 situation. For n=1, the approximation 
is much better, and for the important case where 
n=2, the error introduced is completely negligi- 
ble. 

(1) n=0. The n=0 case corresponds to constant 
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damping. For the two-body interaction Vp'p» we 
consider a separable potential of the form 


V_=-V, le 


b'p - 


pl le, < hw»; 
Voip =0, otherwise. 
With n=0, setting Ap= 0.in the integral equation 
(7) allows us to determine the minimum value of 
qa such that no transition to the superconducting 
state occurs; i.e., T,=0°K. We find a, =hw, 
xexp[-1/N(0)V ], where N(0) is the density of elec- 
tron states of one spin at the Fermi surface. 
This result agrees with that of Suhl.°® : 
(2)n=1. The separable potential of case (1) d 
gives the gap equation 





1 i de 1 
N(O)V . g «(+47 +a)" 















The solution is of the BCS type with Ap=hw)/ 
sinh[(1 + a@?)/N(0)V ] for l€p|<hwo and A, =0 others 
wise. Thus, the inclusion of linear damping modi 
fies the BCS solution only through the introductio 
of a reduced effective coupling constant N(0)V/ ' 
(1+ a). 

(3) n=2. The n=2 case is of interest both for | 
superconductors and He® in which a short-range | 
potential gives rise to damping by exciting parti- | 
cles out of the Fermi sea. For quadratic damp- | 
ing, the T=0°K gap equation is 
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In the case of the superconductor the screened 
Coulomb interaction leads to a@ of the order 1/ SEr, 
indicating that Coulomb damping effects are small 
since matrix elements of V,,, for excitations from 
the Fermi surface of energy SEF are small. 

For He*, damping effects are important in deter- 
mining Ap and the transition temperature to a pos- 
sible superfluid state. Due to the hard-core po- 
tential the integral equation must be solved in co- 
ordinate space. However, we can obtain an esti- 
mate for the reduction of T, due to damping ef- 
fects by considering an effective separable poten- 
tial. The coefficient a can be estimated from the 
thermal conductivity data of Anderson, Salinger, 
and Wheatley.® We find a=9.1x10'* erg™. When 
compared to the He* Fermi temperature of about 
3°K, a=37.6E,, indicating that the imaginary 
and real parts of the excitation energy are equal 
at pg~E,/40. With the separable potential, the 
transition temperature T, is given by 


1 hw, tanh(8 €/2) 
mov), de 7a + 40%) -1)’ 
where 8,.=kT;,. Choosing hwg~Ep, V is deter- 





(10) 


mined by requiring T, in the absence of damping 
to be T,° calculated by Emery and Sessler.’ The 
critical temperature T, in the presence of damp- 
ing is found to be T,=0.0167,°=8x10-*°K. Thus 
it appears that the effects of damping account for 
the absence of a transition to a superfluid state 

in current experiments performed at tempera- 
tures T > 5x10°K, 
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PERSISTENT RING CURRENTS IN AN IDEAL BOSE GAS 


J. M. Blatt 


Bell Telephone Laboratories, Murray Hill, New Jersey, and 
Department of Applied Mathematics, University of New South Wales, Kensington, N. S. W., Australia 
(Received July 10, 1961) 


In their fundamental paper on the theory of elec- 
trical resistivity, Kohn and Luttinger’ considered 
a gas of particles which do not interact with each 
other, but do interact with randomly located scat- 
tering centers. They proved that this system has 
a normal electrical resistivity, and that the sta- 
tistics (Fermi-Dirac or Bose-Einstein) of the par- 
ticles does not make any qualitative difference. 

Their proof breaks down, however, below the 
Bose-Einstein condensation point, where one sin- 
gle-particle quantum state is occupied by a mac- 
roscopic number of particles. We have investi- 
gated this case for a system contained in a ring- 
shaped volume, and find that persistent ring cur- 
rents are set up by the usual method. That is, 
the ring is placed between the poles of a magnet, 
then cooled below the Bose-Einstein transition 
temperature, and finally removed from the ex- 
ternal magnetic field. Within the assumptions 
of the model, the lifetime of these currents is 
infinite, independently of the physical size of 
the ring. 

The phenomenology corresponding to this model 
is similar to, but not identical with, the London? 
equations. In particular, the “flux quantum” hc/e 
enters the equations explicitly. If the “bosons” 
are taken to be electron pairs,® the charge e is 
twice the electronic charge. The relation between 
the flux through the ring, and the “external” 
flux @g, is given in the model by 


&-me, =(@, -me,)/(1+qN,). (1) 


Here ®, is the “flux quantum,” m is the integer 
closest to the ratio 62/1 at the time the system 
is cooled through the transition temperature, N, 
is the number of condensed bosons (electron 
pairs), and qg is the quantity 


q=e" £/(2np)*M, (2) 


where £ is the self-inductance of the ring, p is 
its mean radius, and e and M are the charge and 
mass of the boson. In applications to supercon- 
ductors, e and M should be taken as twice the 
charge and inass of an electron, respectively, 
and in the limit of low temperatures, N, is half 
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the number of conduction electrons in the sample. 

As the ring is cooled down through the transi- 
tion temperature, without being moved, %, re- 
mains constant, but small ring currents are gen- 
erated which make the actual flux ® very close to 
an integral number of flux units, m#,.* When the 
ring is then pulled out of the external magnetic 
field, dg decreases to zero, and Eq. (1) predicts 
a new value for the final flux @ maintained by the 
persistent currents. This final flux may differ 
appreciably from an integral number of flux “quan- 
ta,” the difference being given by 


(@ -m®,)/®, =-m/(1+qN,). (3) 


With the experimental conditions of reference 4, 
m is of order unity and gN, much larger than 1; 
thus quantization should be maintained. But it 
should be possible to construct experimental con- 


ditions so as to observe the predicted deviation (3). 


It is well known’ that the Meissner effect in ac- 
tual superconductors can be understood in precise 
analogy to such a Bose-Einstein gas model, the 
“bosons” being attractively correlated electron 
pairs. It seems likely, therefore, that the model 
calculation of persistent currents reported here 
can also be extended to an appropriate calculation 
for actual superconductors. There are, of course, 
significant differences between this model and ac- 
tual superconductors, the most important differ - 
ence being the complete absence, in the model, 
of an energy gap. We conclude that an energy gap 
is not required for the understanding of persistent 
currents; for if it were, we would not have been 
able to obtain persistent currents in the model. 

The strictly infinite lifetime of ring currents 
in this model can be traced directly to difficulties 
in the foundations of nonequilibrium statistical 
mechanics’; these difficulties are present in the 
Kohn-Luttinger theory, but are concealed in this 
theory by the extensive averaging necessary at 
high temperatures. 

The author is grateful to Dr. P. W. Anderson, 
Dr. L. M. Delves, and Dr. R. Garwin for valuable 
and illuminating discussions, and to Bell Tele- 
phone Laboratories for a visit in the course of 
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which many of these ideas were conceived. 
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MILLIMETER CYCLOTRON RESONANCE EXPERIMENTS IN DIAMOND 


Conrad J. Rauch 
Lincoln Laboratory,* Massachusetts Institute of Technology, Lexington, Massachusetts 
(Received June 29, 1961) 


Cyclotron resonance has been observed in sev- 
eral semiconducting diamonds at 70 kMc/sec and 
at helium temperatures. For magnetic fields be- 
low 33.5 kilo-oersteds, two strong lines have been 
observed having effective masses of (0.70+ 0.01)m, 
and (1.06 + 0.04)m, with little or no anisotropy. The 
experiments will be described in this Letter. 

The 4-mm wave spectrometer system is similar 
to a 2-mm spectrometer described previously,’ 
with the following exceptions. A three-port cir- 
culator? and a commercial 1N2792 crystal detec- 
tor are used to observe the signal reflected from 
the microwave cavity which operates in the TM,,, 
mode and has a loaded @ factor of about 1500. As 
previously, the samples are positioned in the high 
axial electric field of the cavity for cyclotron reso- 
nance. The static magnetic field is obtained from 
a commercial 12-inch electromagnet with a 9/16- 
inch pole face gap. The light modulation of car- 
riers and detection system as described previous- 
ly is used. 

The diamond crystals which were investigated 
are of the IIb variety.* Resonance was observed 
in only a few of the samples, and there appears 
to be a qualitative agreement of crystal perfection 
as determined from x-ray data with the observed 
linewidths in the spectrum. Some samples having 
wT less than unity exhibited the characteristic 
nonresonant microwave magnetoresistance. Fig- 
ure 1 shows a typical experimental trace in one 
of the better samples. 

The curve is a reproduction of a recorder trace 
of the absorption at 69.2 kMc/sec as a function of 
magnetic field. The data were obtained at an ~0.1- 
microwatt rf power level at 1.2°K with a low level 
of light excitation. Lines are observed at 0.70m, 
and at 1.07m, having linewidths given by w1~=13 
and wt=7, respectively. The recorder trace of 
absorption is taken as a linear function of mag- 


net current but the magnetic field is linear only 
through the low-field line. The high-field line 
falsely appears to be asymmetrical and much 
broader than twice the low-field line, but this 
results from a progressive increase in nonline- 
arity of the magnetic field. 

Absorption line intensities decrease upon moving 
the sample to a region of weaker rf electric field, 
thus indicating that cyclotron rather than paramag- 
netic resonance is observed. Detailed anisotropy 
measurements indicate that little or no aniso- 
tropy exists in either resonance within the limits 
of accuracy of the experiment. At 4.2°K, the 
spectrum below 33.5 kilo-oersteds is very similar 
with the exception of unresolved structure on the 
low-field side of the low-mass line. On occasion 
at relatively low rf level, both at 4.2°K and 1.2°K, 
a weak line at about 0.5m, is observed which may 
be a “quantum effect” line. 

An attempt to observe resonance in the conduc- 
tion band by exciting carriers across the gap 
(5.5 ev) with a hydrogen arc has not been suc- 
cessful. 

The carriers which are observed have been ex- 
cited with a tungsten light source which does not 
visibly glow, thus suggesting that the resonances 
are due to carriers in the valence band since the 
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FIG. 1. Experimental trace of cyclotron resonance 
absorption in diamond at 69.2 kMc/sec and 1. 2°K. 
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impurity levels in the IIb’s lie a few tenths of a 
volt above the valence band.* An estimate of the 
spin-orbit splitting A of the order of 0.005 volt 
from atomic data suggests that possibly the con- 
stant-energy surfaces cannot be simply expressed 
by® 


— 4, C2(p 2p 2 2p 2 2p 2))2 
E(k) = Ak? + [B*k*+ C7 hth te tk Me )}, 


which is good for the valence band of germanium® 
and silicon,’ but that higher order terms must be 
taken into account in the secular determinant 

for diamond. It has been suggested® that the low- 
mass line is due to the light hole, and that the 
high-mass line, being isotropic, is a spherical 
energy surface which goes with the split-off band, 
and that the heavy hole should be anisotropic and 
have a mass of ~3m,. Further experiments are 
in progress to clarify the data. 

Thanks are due to Dr. J. F. H. Custers of the 
Diamond Research Laboratory, Johannesburg, 
for providing the variety of samples. R. N. Dex- 
ter of the University of Wisconsin participated in 


a portion of the experiments. The continued in- 
terest and many invaluable discussions with H. J. 
Zeiger are gratefully acknowledged. W. C. Ker- 
nan and F. H. Perry assisted in making the meas- 
urements. 





*Operated with support from the U. S. Army, Navy, 
and Air Force. 
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ELECTRIC -FIELD-INDUCED MODULATION OF THE ABSORPTION DUE TO INTERBAND 
TRANSITIONS OF FREE HOLES IN GERMANIUM 


M. A. C. S. Brown and E. G. S. Paige 


Physics Department, Royal Radar Establishment, Malvern, England 
(Reveived June 2, 1961) 


Provided the structure of the valence band of 
Ge in the vicinity of k=0 is independent of tem- 
perature, then the absorption due to direct inter- 
band transitions’* depends on temperature only 


through the distribution function of the free holes. 


The distribution function is of prime importance 
in determining the absorption, and any agency 
which changes the energy distribution of holes in 
K space will produce a corresponding change in 
absorption. It is possible, therefore, to study 
the effect of a high electric field (E) on the dis- 
tribution function in a direct manner. In prin- 
ciple, if the details of the band structure are 
known, it is possible to relate a photon energy 
with a range (because of warping) of Kk for both 
types of hole and hence to obtain the distribution 
function in k space of both light and heavy holes. 
Such a relationship is shown in Fig. 1 for transi- 
tions to the split-off band, for a band structure 
calculated as described below. The arrows indi- 
cate the value of k* at which the energy of the 
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hole (light or heavy according to band) is 3kT. 

In practice the change in absorption was ob- 
served due to a change in distribution produced 
by the high field; such a differential method is 
inherently more sensitive than a method meas- 
uring the total absorption. In a subsidiary ex- 
periment performed for the purpose of compari- 
son, the distribution function was changed in a 
known way by raising the lattice temperature a 
small amount and the change in absorption ob- 
served (to be referred to as the 5T effect). The 
change in absorption can be calculated from 
Kane’s* theory. Taking the spin-orbit splitting 
as 0.290 ev, and the other band-structure param- 
eters as given by Kane—but modifying the energy 
dependence on k of the split-off band by 30%, as 
Kane suggests, crudely to allow for higher order 
perturbation terms not rigorously included-—the 
change in absorption due to a 13.5% increase of 
temperature of a Maxwellian distribution at 93°K 
(Fig. 1, curve 1) and 293°K (Fig. 2, curve 1) has 
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FIG. 1, Change in absorption with photon energy for 
a lattice temperature of 93°K. Curve 1, theoretical 
6K for transitions between heavy-hole and split-off 
band (6T/T=0.135). Curve 2, field effect for E=170 
vem‘, Curve 3, 5T effect for 5T=12.6°K. Below, 
the relationship between k? and photon energy. 


been calculated for transitions from the heavy- 
hole band to the split-off band. (The free carrier 
concentration was chosen to be 4.4 x10*5 cm™.) 
In the experimental arrangement for measuring 
the field effect, light from a tungsten filament 
lamp was passed through a germanium filter (to 
prevent optical excitation of carriers in the speci- 
men) and focussed on the optically polished ger- 
manium specimen mounted in a low-temperature 
cell. After being passed through a Leiss spec- 
trometer, the light was finally focussed on a 
cooled (77°K) InSb detector. High-field pulses 
of 5-usec duration and of repetition frequency 
1.3 sec“! were applied to the dumbbell-ended 
Specimen. Though it caused a significant rise 
in lattice temperature, a pulse of 5-ysec dura- 
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FIG, 2. Change in absorption with photon energy for 
a lattice temperature of 293°K. Curve 1, theoretical 
6K for transitions between heavy-hole and split-off 
band (6T/T=0.135). Curve 2, 6K x5 for field effect 
at E=1480 vem™, Curve 3, 5K x3 for 6T effect, 
6T =13°K. 


tion was necessary to enable a saturation of the 
field effect to be observed since the time con- 
stant of the detector was 1 psec. However, it 
was always possible to distinguish and eliminate 
thermal effects because of the relatively long 
thermal time constant of the specimen. Deter- 
mination of the total intensity of the light trans- 
mitted enable the change in absorption (5K) to 
be calculated. 

The subsidiary experiment was performed using 
the same optical arrangement. The change in lat- 
tice temperature was produced by passing a di- 
rect current through the specimen for 2 sec.. A 
steady state was reached during this period, the 
temperature increase being indicated by a thermo- 
couple of suitable response time and in good ther- 
mal contact with the specimen. The beam was 
stopped down so that only a volume of the crystal 
close to the thermocouple was “active.” 

Measurements have been made at various lat- 
tice temperatures between 90°K and 320°K on a 
p-type Ge sample containing 4.4 x10'* cm* free 
carriers. Curve 2 of Fig. 1 shows 6K for a field 
of 170 v cm™ and curve 3 shows the 5T effect for 
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5T/T=0.135, both measured with the lattice at 
93°K. A similar pair of results are shown in Fig. 
2 for a lattice temperature of 293°K. The change 
of absorption due to transitions of heavy holes to 
the split-off band can readily be identified by com- 
parison with the calculated curve; a maximum at 
about 0.27 ev is associated with transition of the 
light holes. 

Some of the salient points which emerge from 
these results and from the velocity field curves 
for the same sample are: 

(1) At a lattice temperature of 93°K, a com- 
parison between the field effect and the 57 effect 
shows that in electric fields up to 350 v cm™ the 
distribution remains Maxwellian, within the ac- 
curacy of the measurement. Figure 3 shows the 
relationship between 5T and E obtained from the 
comparison. The calculated change in absorption 
of heavy holes is in reasonable agreement with 
the observed values both with regard to shape and 
magnitude. 

(2) At a lattice temperature of 293°K there is 
a marked difference between the field effect, the 
5T effect, and the computed 6K. If the low-tem- 
perature band-structure parameters are correct, 
then it is possible that the discrepancy between 
the predicted curve and the 5T effect may be due 
(a) to changes of the band structure with temper- 
ature to which the calculated and observed results 
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FIG. 3. The increment of hole temperature with 
electric field strength for a lattice temperature of 
93°K, Maximum error in 6T is +15 %. 
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are sensitive, or (b) to indirect transitions. Al- 
ternative (a) is favored since for indirect transi- 
tions 6Kjng/5T is always positive, and 5Kjng has 
been estimated to be less than 20% of the value 
for direct transitions. In an electric field the de- 
parture of the distribution from Maxwellian at 
this higher lattice temperature can be understood 
since the cross section for hole-hole scattering 
falls* while the cross section for phonon scatter- 
ing rises with increasing temperature. 

(3) For given scattering mechanisms, Fig. 3 
enables the field dependence of the mobility to be 
predicted and, from the energy balance equation, 
a value of (=jop/=lac)* can be determined if the 
same value is assumed for all transitions. Taking 
Brown and Bray’s® values of the parameters, 
which occur in the mobility formula, we find (i) 
that the calculated ratio of zero-field mobilities 
at 93°K and 293°K agrees with the observed ratio 
to within 4%; (ii) that a dependence of lattice mo- 
bility, 4,, on the hole temperature, T;, of uz 
«Tp -6, §~1.3,"has to be assumed to account for 
the observed dependence of mobility on field; and 
(iii) that the mean value of (=jop/=lac)” is 1.6+ 0.3, 
which is approximately half the value required if 
optical and acoustical phonon scattering are re- 
sponsible for the .; « T~*** law, and implies that 
the normal form of optical and acoustical phonon 
scattering is not entirely responsible for the high 
value of 6. (Calculations we have made show that 
the nonparabolic shape of the light-hole band® 
makes an insignificant contribution to the temper- 
ature dependence of 1; .) 

In conclusion we should like to point out that the 
theoretical response time of the modulation is the 
energy relaxation time 7, of the free holes. At 
room temperature 7, ~10™ sec,® and is not more 
than 107*° sec at 93°K. Since substantial depths 
of modulation (>10%) can be achieved with the 
lattice at 93°K, the field effect we observe could 
form the basis of a modulator capable of modu- 
lating a light beam at very high frequencies. 
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HOT ELECTRON EMISSION FROM SILICON pn JUNCTIONS PARALLEL TO THE SURFACE* 


J. L. Moll, N. I. Meyer,t and D. J. Bartelink 
Solid-State Electronics Laboratories, Stanford University, Stanford, California 
(Received June 19, 1961) 


A number of papers have been published on 
emission of electrons from semiconductor pn 
junctions into vacuum.** It is a common feature 
of these works that the emitted electrons have 
come from regions where the junction intersected 
the surface. This emission geometry limits the 
magnitude of the emitted current and makes a 
quantitative interpretation of the experimental 
results difficult. An attempt to observe emis- 
sion from a pn junction parallel to the emitting 
surface was mentioned by one of the authors.* 
The attempt did not succeed—apparently because 
the neutral n layer between the junction and the 
emitting surface was too thick (~ 1400 A). 

In this Letter we present some results on elec- 
tron emission from silicon pm junctions parallel 
to the emitting surface. The junctions are made 
from cubes of 0.02 ohm cm p-type silicon approx- 
imately one millimeter on the edge. Four faces 
were diffused with thick n+ regions whose purpose 
was to make contact to one face on which a thin 
n-type region was diffused. The thin layer was 
about 1000 A thick. One face was used to con- 
tact the p-type body, thus allowing bias to be 
applied to the junction. Typical breakdown volt- 
ages were 10 volts for the thick junction and 4.5 
to 6 volts for the thin junction. 

The emission current was measured in a de- 
mountable vacuum system with the device posi- 
tioned at the center of a hemispherical collector. 
A small hole was drilled in the collector to allow 
light to be shone on the emitting surface. At pre- 
breakdown reverse bias voltages the junction cur- 
rent was typically raised from about 50 pa to 1 
ma by the light. The emission current due to 
these photoelectrons is called photoenhanced 
emission current. Without the light, the junction 
voltage had to be raised to the breakdown region 
in order to give rise to a measurable emission 
current with the sensitivity available (~107** 
amp). Prior to each measurement the junctions 
were etched in HF acid to remove silica and 
placed in the vacuum system without exposure 
to air. 

A method was developed for thinning down the 
n region in well-defined steps. It consists of 
growing layers of SiO, on the crystal in boiling 
water. The SiO, layer is subsequently removed 
by etching in HF. The rate of removal was cal- 


ibrated on uniformly doped silicon crystals by a 
weighing technique and was found to attain a value 
of about 34 A per boiling. This is a convenient 
value with a view to constructing the emission 
devices. When the boiling water technique was 
applied to crystals with a pn junction close to 

the surface, it was found that the rate of removal 
was reduced to about 5 A per boiling when the 
surface came to within about 200 A of the field 
region of the pm junction. 

The emission current is shown in Fig. l asa 
function of the number of boiling water, HF steps 
and is seen to increase exponentially with amount 
of material removed. The exponential increase 
levels off between 18 to 25 steps of boiling. This 
is interpreted as being due to the reduction in the 
rate of removal as discussed above. 

Due to the relatively simple geometry it has 
been possible to set up a mathematical model of 
the device and to derive theoretical expressions 
for the emitted electron current. The theoretical 
results are based on an approximate solution of 
Boltzmann’s transport equation taking into ac- 
count the important types of interaction between 
the hot electrons and the crystal lattice.’ These 
interactions are scattering by optical phonons 
and secondary ionization. The theoretical ex- 
pressions include some basic parameters de- 
scribing these interaction processes. The values 
of the interaction parameters are presently sub- 
ject to controversial discussions.* Analysis of 
our experimental results seems to resolve some 
of the questions concerning the appropriate val- 
ues to be used for the interaction parameters. 

The energy distribution of the hot electrons 
after they have been accelerated in the strong 
field of the pn junction and subsequently diffused 
to the emitting surface through the neutral region 
is mainly determined by the following quantities: 

The field E in the pn junction, the width w of 
the pn junction, the thickness L of the neutral 
region, the mean free path 1, for optical phonon 
scattering, the mean free path /; for ionization, 
the energy 6 , of the optical phonons, and the 
threshold energy &; for ionization. The param- 
eters E, w, and L may be varied independently 
(within certain limits) by choosing the appro- 
priate device parameters (doping levels, im- 
purity profile, etc.). The quantities l,, lj, &,, 
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and &; are characteristic of the material in ques- 
tion. For silicon &,=0.063 ev as found by neu- 
tron diffraction experiments.® The values of the 
remaining three parameters are subject to dis- 
cussion as mentioned above. 

The magnitude of the emitted current depends 
in addition to the above parameters on the elec- 
tron affinity §.,. The value of 6.» has been de- 
termined from the retarding potential measure- 
ments for an untreated silicon surface. 





where 
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The values of &;, 1,, and 1; enter into the 
theory in a manner such that one single meas- 
urement can be fitted by a wide range of values. 
However, by considering retarding field meas- 
urements, variation of current with thickness, 
and total current, it is possible to narrow down 
the possible range of values. 

The theoretical treatment gives the following 
dependence’ of emission current on L: 


‘.* exp(-L/a), 


(1) 
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(2) 
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and T,1 is an equivalent temperature describing 
the energy distribution of the electrons with en- 


ergies less than 6, at the end of the field region. 


The values of A found from Fig. 1 are A,=20A 
and A,=35 A. For the two extreme theories dis- 
cussed by Shockley,® we find the following limits 
for 1; by using the above experimental values for 
, together with the theoretical expression (2): 


35 A<1.<60A for 1.«l ; (3) 
i . 
300 A</.<500 A for 1.=181 , 
i i Y 
hence 16 aie 28 A. (4) 
The first case corresponds to the results of 


Wolff?® and others, while the second case cor- 
responds to Shockley’s theory.® 


Figure 2 shows the results of a retarding po- 
tential experiment. From the value of the cutoff 
voltage Vem, the electron affinity has been de- 
termined from 


é. =3.7 ev, (5) 


Pi aad 


appl ; @V om 
since the built-in potential is approximately equal 
to the energy gap & g for the junctions in ques- 
tion. 

It is found theoretically that the high-energy 
part of the distribution function is approximately 
Maxwellian with an effective temperature: 


L +l, ‘. kT 
Tag? y/ r Bites < ] 
i yr t 


The value found from Fig. 2 is Tg9 =0.5 volt. 


- (6) 
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FIG, 2. Photoenhanced emis- 
sion current as a function of en- 
ergy. J(&) is the contribution 
to the emitted current from 
electrons with energies between 
& and &+d&(measured in ar- 
bitrary units). The result is 
obtained by differentiating an 
experimental curve for the mag- 
nitude of the total emitted cur- 
rent as a function of the applied 
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Using the values from (3) together with (6) gives 
rise to a value of T,9 which is an order of mag- 
nitude too small. For the values in (4), To2 
comes out about a factor of two too small. 

Finally, the experimental value of the total 
emitted current was compared to the theoretical 
results. Using the values in (3), the calculated 
emission current was found to be several orders 
of magnitude smaller than the measured current. 
The same applies to the theoretical emission 
current corresponding to the lower limit in (4) 
while the upper limit gave reasonable agreement 
between calculated and measured results. If we 
assume a smaller ratio of 1;/l, than 18, all of 
the data can be fitted with 7;~100 A and 1,=50 A. 
One of these parameters can be changed by + 50% 
and a fit still obtained by adjusting the other one. 
As I; is increased, Ll, must be decreased to re- 
tain the fit. 

It is concluded that the low values of 1; postu- 
lated by Wolff and others are not consistent with 
our results. The actual result does not satisfy 
either Shockley’s assumption that 1; >1, or 
Wolff’s assumption that 1, >1;, but seems to 
favor Shockley’s picture. There is an inherent 
experimental error in the retarding potential 
measurement in that the potential across the 
emitting face may vary by a few tenths of a volt 
giving too high a value for Tz9. If Tgo is reduced 
by a factor of 2, Shockley’s picture fits fairly 
well but with a surprisingly small mean free path 


for optical phonon excitation. It would be of great 
interest to study the electrons which have lower 
energies than the 3.7 ev corresponding to the 
electron affinity of the untreated silicon surface. 
This may be accomplished, e.g., by coating the 
surface by Cs and thereby lowering Ss Prelim- 
inary results on devices of the type described 
above with a Cs layer on the surface have been 
obtained. The value of &,, was about 1.8 ev and 
the maximum emission current was 0.02 wa. No 
detailed comparison with theory has been made 
so far, as it is believed that the magnitude of the 
emission current is influenced by the presence 

of an oxide layer between the silicon and the Cs 
layer. 
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LINEAR EFFECT OF APPLIED ELECTRIC FIELD ON MAGNETIC HYPERFINE INTERACTION 


N. Bloembergen 


Harvard University, Cambridge, Massachusetts 
(Received July 10, 1961) 


It has been shown that an electric field will 
change the quadrupole coupling of a nucleus 
which is not at a center of inversion symmetry.** 
The purpose of this note is to show that the mag- 
netic hyperfine interactions, either isotropic or 
anisotropic, can also be linear functions of the 
applied electric field. Consider first the isotropic 
or contact interaction AI-S. The variation in the 
scalar quantity A may be written as AA=F-E. 
There are only a few physical examples of such a 
linear relationship between a polar vector and a 





(016 - F,) la) E15, 10) 


scalar. The best known is the pyroelectric effect. 
The properties of F in different point group sym- 
metries have been described by Voigt.* In our 
case the symmetry at the nuclear site should be 
taken. The relationship between the anisotropic 
hyperfine tensor and the electric field is described 
by a third order tensor. The phenomenology of 
this case has been discussed before. 

A quite general expression for the effect which 
is bilinear in the contact hyperfine interaction 
and the electric field can be written down by 
means of second order perturbation theory: 


+complex conj. (1) 
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The electric field operator Hg = -)-pepExp is odd 
and 6(¢ - Fy) is even with respect to the nucleus 

as origin, Ty =0. A nonvanishing result is ob- 
tained only if the wave functions of both the ground 
state ]0) and the excited state |i) contain some 

s character, and in addition are hybridized by 
some p-wave functions. A necessary condition 

is that the nucleus is not at a center of inversion 
symmetry. 

The effects described above are observable, 
for example, in MnF,. The F’® resonance has 
been thoroughly investigated in both the paramag- 
netic and antiferromagnetic regions.®*® The te- 
tragonal crystal structure is well known and has 
the symmetry Dqg,. All Mn** ions are at a center 
of symmetry. There will be no linear effect of 
an electric field on the dimension of the unit cell 
nor on the magnetic properties of the Mn** sub- 
lattices. The symmetry at the F’® nuclei at posi- 
tions +(u,u,0), +(u, -u, 4), is, however, ortho- 
rhombic, mm2 or Co,. The polar vector F has 
one nonvanishing component. It has the same 
magnitude for each of the four F~ ions in the unit 
cell, but points in the four directions [110], [110], 
[110], and [110]. 

Equation (1) is not suitable for explicit calcula- 
tions, because the wave functions of the excited 
states are not known. The magnitude of the effect 
is determined, as in the case of quadrupole inter - 
action, by the change in interionic distances and 
by the change in covalent character and s-p hy- 
bridization on application of an electric field. 

The change in interionic distances gives a large 
effect in the mostly ionic crystal of MnF,. Bene- 
dek and Kushida® have determined the variation 
of the contact interaction A with hydrostatic pres- 
sure. In their case the dimensions of the unit 
cell, a and c, are varied, whereas u was essen- 
tially constant. If an electric field is applied, u 
changes, but a and c are constant. The variation 
in u can be estimated from the contribution of the 
ionic displacements to the dielectric polarization, 


e mS 
Aary= Auav2 = rhe E. (2) 
The ad hoc assumption has been made that the 
force constants on each F” ion are isotropic in 
the (001) plane. €,, the dielectric constant, and 
n,, the index of refraction, are assumed to be 
equal to the values in CaF,, € =6.8 and n=1.44. 
N is the number of F~ ions/cm*; ry is the varia- 
tion in distance to the nearest Mn** ion. The 
distance to the two other nearest neighbor Mn** 
ions of type I varies by Ary= Ary sind =0.63Arq7. 
Take the Mn** ion at (0, 0, 0) with spin down. 
Consider an electric field of 10* volts/cm applied 
along [110]. The F’® nucleus at (u, u, 0) will move 
closer to the Mn** with spin down and away from 
two Mn** ions with spin up. The local magnetic 
fields and their variation due to displacement 
are listed in the first two columns of Table I, 
for complete spin alignment on the antiferromag- 
netic sublattices. The variations for the F’® nu- 
cleus at (-u, -u, 0) have the opposite sign. The 
experimental results of Benedek and Kushida have 
been used to calculate (8H /au). It is not necessary 
to make theoretical assumptions about the nature 
of the bond at this stage. Either the exchange cor- 
relation in overlapping ionic orbitals® or an admix- 
ture of covalent orbitals®’’ will give this variation 
of the contact interaction with interionic distance. 
In covalent orbitals there is, however, another 
contribution given in the last column of Table I. 
The electric field will change the degree of co- 
valency A” and the degree of sp hybridization 
s*=1-p?, at constant interatomic distance.® 
Although complete expressions are obtainable 
with the Ritz variational method, the following 
approximations will suffice for our present pur- 


pose: 





(i - “(x,t -s*)¥? 


aE) #£4W..-W 
u 
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2p 862s 


Table I. Contributions to the magnetic field acting on the F'® nucleus at (u,u, 0) in antiferromagnetic MnF, with 
complete spin alignment, and their variation in an electric field along [110]; E = 10‘ volts/cm. 








A, (0H, /8u)(8u/dE)E (a4, /2E) E 
(kilogauss) (gauss) (gauss) 
Dipolar 12.50 +0.425 0 
Contact hyperfine 25.7 -13.6 -14 
Anisotropic hyperfine 1.66 eee cee 
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The factor Ejo-/E =3(€ +2) =2.9 relates the local 
field at the position of the bond to the applied 
field E. The Lorentz relationship has been as- 
sumed, although this, of course, is not strictly 
valid at a noncubic site. The energy difference 
between the 2s and 2p level of the fluorine atom 
is Wop -W9,~20 ev; eX 5p) is the dipole mo- 
ment matrix element connecting these states, 
which may be calculated with Slater-type orbitals. 
Wnn -WF is the energy necessary to transfer an 
electron from a F~ ion to the Mn** ion in the lat- 
tice, and is rather arbitrarily taken as 5 ev. The 
contact hyperfine interaction is proportional to 
2A;*s;’, where the summation is over the three 
nearest Mn neighbors. Note that the degree of 
covalency and s character increases if the F 
atom moves toward the Mn atom. The total rel- 
ative variation of the hyperfine interaction at 
constant u is approximately 


(AH /H) | =[2(A4A/d) + 2As/s](1 +2 sind). (5) 


Take s=4 corresponding to 25% s character, 
while 4?=0.02. It turns out that the variation 

in A contributes an order of magnitude more 
than the variation in s. Combination of Eqs. 
(3)-(5) gives the value -14 gauss in the last col- 
umn. 

One finds from Table I that the total relative 
variation of the F*® resonance frequency in anti- 
ferromagnetic MnF, is |Avg|/vg=0.67x10~ for 
E =10* volts/cm along [110]. This corresponds 
to Avg = 108 kc/sec since vg=160 Mc/sec at 
4.2°K. Since the linewidth is only 50 kc/sec, 
the effect should be readily observable. Experi- 
ments are at present under way. They should 
give valuable information about the covalent con- 
tribution and test the correctness of the theoreti- 
cal assumptions. 

The p-type or anisotropic hyperfine interaction 
depends sensitively not only on the o orbital con- 
sidered above, but also on the 7 and 7’ orbitals.” 
A calculation of its dependence on electric field 
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will be postponed until a more complete account 
can be given. The contribution to the F’® reso- 
nance frequency is estimated to be at least one 
order of magnitude smaller than the other entries 
in Table I. 

The dipolar and anisotropic hyperfine interac- 
tion also produce a magnetic field H, in the [110] 
direction at right angles to the z axis of quantiza- 
tion of the F*® spin. Consider the dipolar term, 


“88 y HI = -§sindcos6 
x r*g8g,8,(S_)" +I). (6) 


One finds AH, = -0.6 gauss if an electric field of 
10* volts/cm is applied at 4.2°K. If the electric 
field is alternating at the resonant frequency of 
the F’® spins, it will induce nuclear spin transi- 
tions. At 4.2°K the F’® resonance can be satura- 
ted electrically. This situation is rather unusual, 
in that a spin /=4 is “tickled” by an applied elec- 
tric field. 

Finally it is pointed out that the isomer shift,® 
which is also proportional to |¥?(0)|, will be a 
function of the electric field, if the nucleus is at 
a suitable lattice site. Since the relative varia- 
tion in |¥7(0)| of the valence electrons is about 
1:10° for 10° volts/cm, the effect will be near 
the limit of detectability for the Missbauer ef- 
fect. 
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IMAGE OF THE FERMI SURFACE IN THE LATTICE VIBRATIONS OF LEAD 


B. N. Brockhouse, K. R. Rao,* and A. D. B. Woods 


Neutron Physics Branch, Atomic Energy of Canada Limited, Chalk River, Ontario, Canada 
(Received June 26, 1961) 


On theoretical grounds Kohn’ has proposed that 
there should exist positions in reciprocal space 
at which the slopes of the dispersion curves of 
the lattice vibrations in metals become infinite. 
The positions are related to the Fermi surface 
of the electrons, and the infinite slope to the 
sharpness of the Fermi surface. These Kohn 
anomalies have not previously been observed 
experimentally for certain, and recent calcula- 
tions? suggest that they may be quite small, at 
least for some materials. We believe that we 
have observed the anomalies in a single cyrstal 
of lead. 

Previous experiments’® on lead at this laboratory 
yielded dispersion curves which could be inter- 
preted on the Born-von Karman theory of lattice 
dynamics only if it were assumed that there exist 
long-range forces between ions in the crystal, the 
forces often varying in sign. These results were 
consistent with the existence of Kohn anomalies, 
but direct evidence was weak. 

The new experiments were carried out using the 
triple-axis crystal spectrometer in such a way 
that, for each energy distribution, the momentum 
transfer, hQ=h(k, -k’), remained fixed.* The cen- 
ters of the neutron groups observed give the fre- 
quencies v of the lattice vibrations (of wave vector 


Q in the extended zone scheme) through the rela- 
tion hv=E,-E’. Inthe expressions k, and k’ and 
E, and E’ are the incoming and outgoing neutron 
wave vectors and energies, respectively. 

Figure 1(a) shows a series of neutron groups 
obtained as the terminus of Q moves along the 
line Q=(21/a)[t, £, £ ] between the reciprocal lat- 
tice points (1,1, 1) and (2, 2,2) in steps of 0.025 
in ¢. (a is the lattice constant.) The neutron 
groups observed correspond to longitudinal pho- 
nons whose reduced wave vector q is in the [111] 
direction of the Brillouin zone. The centers of 
the groups are shown plotted in Fig. 1(b) as a 
function of (a/2n)Q. Near the zone boundary a 
sharp anomaly occurs on both sides of the zone 
boundary, as it should, and three independent 
series of measurements at different energies 
gave almost identical results. 

For a spherical Fermi surface the positions of 
the anomalies are given’ by the equation 


2K re |2n7+ql, 
where 7 is a vector of the reciprocal lattice and 
KF is the radius of the Fermi surface. In Fig. 2 
the (110) plane of reciprocal space is shown, each 


reciprocal lattice point, a7, surrounded by its 
zone. A spherical Fermi surface containing four 
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FIG. 2 The (110) plane of the reciprocal lattice of 
(face-centered cubic) lead, each point surrounded by 
its zone. A spherical Fermi surface containing four 
free electrons (extended zone scheme) intersects the 
plane as shown by the circle about the point (0,0,0). 
In the reduced zone scheme this spherical surface gives 
rise to the surface shown about (1,1,1) in the second 
zone, and to the surfaces shown about (1,1,1) in the 
third zone. Possible effects of gap formation on the 
extended zone surface are indicated by bold lines. The 
Kohn construction for free electrons is shown by the 
circular arcs AECL FIDG and CJDBKH. The points A, 
B---L indicate positions of Kohn anomalies for free elec- 
trons in the symmetric directions [100], [110], and 
{111}, shown by dashed lines. 


electrons®*® is drawn around the point az = (0, 0, 0). 
In the reduced zone scheme this Fermi surface 
gives a full first zone, a second zone with a pock- 
et of holes as shown around the point (1, 1,1), and 
a complicated structure in the third zone shown 

in cross section around the point (1,1, 1). The 
locus of the Kohn anomalies connected with the 
point (0, 0,0) is shown as the circular arc 
AECLFIDG. Other types of anomalies also occur, 
as indicated by the arc CJDBKH which is connec- 
ted with the point (2, 0, 0) [a point not in the (1, I, 0) 
plane]. The positions designated by the letters 

A, B, +++, L are those for which anomalies in 
the symmetric [100], [110], and [111] direc- 
tions of the reduced zone occur for a spherical 
Fermi surface. Some of the anomalies have a 
multiplicity greater than unity, because they can 
be produced from several reciprocal lattice vec- 
tors (e.g., E with multiplicity 3, C with 4, H with 
2, etc.). Some of the anomalies can occur only 

in longitudinal modes (e.g., A, F, G); these have 
unit multiplicity. The signs of the anomalies (i.e., 
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are, 4 =+e or -%) can be assigned by the plausible 
rule that the frequency increases anomalously 
when q passes from the interior to the exterior 

of the Kohn sphere since the electron gas then 
“stiffens.” Thus A, E, F, and / have positive 
signs, while C, D, and G have negative signs. 

The spherical Fermi surface applies to ideally 
free electrons. In fact energy gaps will be pro- 
duced at Bragg reflection planes, as illustrated 
by the modifications to the spherical Fermi sur - 
face shown as bold lines in Fig. 2. These gaps 
round off the Fermi surfaces in the reduced 
zone, and render the Kohn construction’ more 
complicated. In particular, only points A, F, G, 
and possibly J, are certain to be unchanged in a 
qualitative sense by the passage to the real Fermi 
surface. 

It will be observed that point F, expected to oc- 
cur only in the longitudinal mode, corresponds 
very closely to the anomaly in Fig. 1(b). However, 
this assignment is ambiguous because of possible 
interference with point E. 

In Fig. 3 the longitudinal branch for the [110] 
direction is shown. These measurements were 
made along three separate lines through recipro- 
cal space; from (2, 2,0) to (1,1, 0), from (2, 2, 0) 
to (3,3, 0), and from (1,1, 1) to (2,2,1). Two dif- 
ferent monochromator planes were used from 
(2, 2,0) to (1,1,0). An anomaly, which can be 
identified with point G, appears at aqg/27= 0.40. 
The position of the anomaly yields the result that 
the radius of the Fermi surface in the [110] direc- 
tion (extended zone scheme) is about 1.21(27/a), 
instead of the value 1.24(27/a) holding for free 
electrons. Thus the width of the arms in the third 
zone (reduced zone scheme) is smaller by about 
15% than for free electrons, in qualitative agree- 
ment with deductions by Gold’ from his measure- 
ments on the de Haas—van Alphen effect. [Note 
added: From inspection of Fig. 3, Kohn and Woll 
(private communication) would place the anomaly 
at ~ 0.47, and thus make the width about one third 
smaller than for free electrons. | 

Many other features of the unusual dispersion 
curves of lead® may be interpreted in terms of 
the Kohn effect. Of particular interest is the 
spectacular drop in frequency found for the longi- 
tudinal branch in the [100] direction near the zone 
boundary, which may be connected with anomaly 
C. This anomaly has fourfold multiplicity and 
therefore should be exceptionally strong. However, 
details are not clear since this point is strongly 
influenced by gap formation at the (002) Bragg 
plane, as indicated in Fig. 2. 
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| FIG, 3. The longitudinal 
branch in the [110] direction of 

the reduced zone. The branch 
is continued beyond the zone 
‘ 7 boundary to the point (1,1,0). 
° The continuation, lying on the 

™ 7 square face of the Brillouin 
zone, has the same symmetry 
as the [110] direction within the 
x zone. The anomaly assigned as 
G is indicated by the solid verti- 
cal arrow at ga/2n~0.40. For 
: free electrons it would appear 
at 0.35 as indicated by the 
dashed arrow. Other possible 
anomalies can also be seen. 
Measurements in a single series 
are thought to have relative er- 
rors about the size of the points. 
The scatter in points belonging 
to different series may be indic- 
ative of the absolute errors. 
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As mentioned earlier, recent calculations? have 
indicated that the anomalies should be quite small. 
The strength of the effect in lead is probably con- 
nected with the strong phonon-electron interaction 
inthis material, as indicated for example by the 
comparatively high superconducting transition 
temperature. This is supported by the fact that 
the Kohn effect has not been detected in similar 
measurements on a single crystal of sodium.’ 

The anomalies are, of course, smeared out by 
resolution in 4 space. (Indeed, Fermi surfaces 
may be more easily delineated by means of x-ray 
intensity measurements than by neutron measure- 
ments, because of the better resolution in 4 space 
available with x rays.) However, even under good 
resolution the anomalies will appear to be smeared 
because of the logarithmic character of the infinity 
in the slope, which makes exceedingly small the 
range in vy over which |grad=v| is large. 

Dr. G. Caglioti participated in an earlier, incon- 
Clusive, attempt to detect the effect in lead. E. A. 
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Glaser gave valuable technical assistance. 
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DISSOCIATION OF H,” IONS BY A MAGNETIC FIELD* 


Selig Kaplan, George A. Paulikas, and Robert V. Pyle 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received July 10, 1961) 


The dissociation of molecular ions by an ex- 
ternally applied electric field, E, or an equiva- 
lent electric field, E’=v xB, due to motion 
through a magnetic field has been calculated by 
Hiskes.' Qualitatively, it is predicted that the 
external fields warp the nuclear potential in such 
a way that the higher vibrational states become 
unstable and dissociation results. This effect 
has been observed recently by Riviere and Sweet- 
man,? who dissociated H,* and H,* ions of approxi- 
mately 2-Mev energy with electric fields of up to 
5x10° v/cm. The measurements presented here 
demonstrate the breakup of molecular hydrogen 
ions of 10 Mev per nucleon by static magnetic 
fields. 

The experimental arrangement is shown in 
Fig. 1. Ions produced in a P.I.G. discharge were 
accelerated in the Berkeley heavy-ion linear ac- 
celerator (Hilac). H,* and H,* ions were continu- 
ously accelerated to full energy, and, alternative- 
ly, full-energy H,* ions were obtained by partially 
stripping 3-Mev H,* ions in a jet of Hg vapor and 
then accelerating the resulting H,* ions to 20 Mev. 
The reason for producing Ny by the latter tech- 
nique is that such ions might be left in more high- 
ly excited states than those which came directly 
from the ion source, as demonstrated by Riviere 
and Sweetman. The ions emerged from the Hilac 
accompanied by other particles formed by break- 
up on the residual gas in the vacuum system. The 
beam was collimated to % x % in., and the various 
particles were then separated by bending the use- 
ful beam by 15 deg in a dc magnetic field. No 
further collimation was required, thus eliminat- 
ing background from breakup at slit edges. 
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FIG, 1. Experimental arrangement, 
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The dissociation and analysis were obtained in 
a dec magnet with a 4-in. gap. The gap in the first 
1.5 in. of the magnet could be reduced to 1 in., 
thereby producing a maximum magnet field of 
about 19 kilogauss. The beam required about 10~° 
second to traverse this region. 

The results reported here were obtained by acti- 
vating 0.010-in. thick Cu foils, placed in the ana- 
lyzing magnet perpendicular to the direction of 
the incident beam. The foil holder was built to 
serve as a Faraday cup, and was used in the beam- 
monitoring system. The spatial distribution and 
the relative intensities of the beam particles (i.e., 
H°, H,*, H,*, etc.) were determined by counting 
the Zn® activity (7,,.=38 min) produced in the foil. 

Briefly, the scanning was done as follows: The 
bombarded Cu foil was wrapped around a rotating 
wheel whose instantaneous position was synchro- 
nized with the voltage of a linear ramp. A Nal 
crystal viewed the wheel through a slit in a W and 
Pb shield. Whenever a 0.51-Mev annihilation y 
ray from Zn® decay was detected in the crystal, 
a pulse of height equal to the ramp voltage was 
generated. These pulses were put into a Penco 
100-channel pulse-height analyzer, and a presen- 
tation of count rate vs position on the foil was ob- 
tained. 

The dissociation caused by the residual gas in 
the vacuum system at our base pressure of about 
6x107° mm Hg was determined by raising the 
pressure in several steps by adding N,, meas- 
uring the dissociation, and extrapolating to the 
base pressure. 

The H,* ions obtained by stripping H,* were 
partially dissociated by the magnetic field, with 
a rapid rise at about 12 kilogauss. This magnetic 
field corresponds to approximately 510° v/cm 
(Fig. 2). Our result is consistent with the disso- 
ciation of the v=17 vibrational state as calculated 
by Hiskes for nonrotating molecules. The verti- 
cal bands marked 107° sec and 107 sec indicate 
the calculated fields necessary for barrier pene- 
tration in these times, with the widths of the bands 
corresponding to the uncertainties. The field in 
the 15-deg bending magnet was sufficiently high 
to have already dissociated the v = 18 vibrational 
state if it were populated. Conversely, the maxi- 
mum stripping field was somewhat too low to af- 
fect the v=16 state. 
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FIG. 2. Fraction of the beam dissociated by magnetic 
fields (present experiment) and electric fields (Riviere 
and Sweetman). a, magnetic, H,* from H;*; @, mag- 
netic, H,* direct; A, electric, H,* from H,*; 0, elec- 
tric, H,“ direct. A background-gas breakup correspond- 
ing to a fractional dissociation of (0.26 +0. 04) x107° 
has been subtracted from the magnetic dissociation 
data. The standard deviations shown include counting 
statistics, an estimated uncertainty of 15 % in the ex- 
trapolated dissociation by the residual gas in the vac- 
wm system, and an uncertainty of 5% in the calibra- 
tions of foil activity vs ion current. 


No indication of magnetic field dissociation was 
obtained for the directly accelerated H,* and H,* 
ions. We have no theoretical information concern- 
ing the H,* ion. The absence of dissociation of 
directly accelerated H,* ions is assumed to be due 
to the method of formation and to possible de-exci- 
tation in the ion source. 

With more suitable choices of the kind of ion, 
energy, and magnetic field strength, dissociation 
by a magnetic field might provide a mechanism 
for trapping ions in a highly evacuated controlled 
thermonuclear device, and provide a diagnostic 
tool for the study of molecular ion formation and 
excitation. We are extending this investigation 
to other molecular ions and to neutral atoms in 
similar and higher magnetic fields. 

We wish to express our appreciation to Dr. C. M. 
Van Atta for supporting and encouraging this re- 
search, J. Warren Stearns and Henry F. Rugge 
for helping with the measurements, Dr. Edward 
L. Hubbard, Duane A. Spence, and other members 
of the Hilac staff for their assistance, John R. 
Meneghetti for much of the mechanical construc- 
tion, and Dr. John R. Hiskes for many interesting 
and useful theoretical discussions. One of the au- 
thors (S. K.) would like to thank Professor Burton 
J. Moyer for the support and interest that enabled 
him to participate in this work. 
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ROTATIONAL BANDS IN Ne” 


A. E. Litherland, J. A. Kuehner, H. E. Gove, M. A. Clark, and E. Almgqvist 
Atomic Energy of Canada Limited, Chalk River, Ontario, Canada 
(Received June 30, 1961) 


In recent years evidence has been obtained»? 
which strongly suggests that the low-lying levels 
of Mg”*, Mg?®, and Al*® can be arranged into rota- 
tion-like bands. In this Letter we discuss new 
evidence which suggests that the low-lying levels 
of Ne*° can also be arranged into overlapping rota- 
tional bands. A separation of the level spectrum 
of Ne*° into rotational bands is shown in Fig. 1. 
Ne” is of great interest in this respect because 
it is close enough to the closed shell at O'* for 
detailed shell-model calculations to be feasible.*® 

Last year Broude and Gove?* showed that the 
spins of the 4.25- and 4.97-Mev levels in Ne”° 
were 4 and 2 and that the quadrupole to dipole 
amplitude ratio for the cascade transition from 
the 4.97- to the 1.63-Mev level was +0.08. These 
results have been confirmed recently and the pari- 
ties of the levels measured.°~’ The parity of the 
4.97-Mev level was found unexpectedly to be nega- 
tive and the astrophysical implications of this 
measurement have been discussed recently.’ 
Measurements on the lifetimes of the levels in 
Ne”° using the Doppler shift attenuation method®,® 
and Coulomb excitation’® have indicated that the 
electric quadrupole transitions from the 1.63-Mev 
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FIG, 1. A separation of the level spectrum of Ne” 
into rotational bands. A level in Ne” above about 4.75 
Mev can break up into O"* and an alpha particle if its 
parity is (-1)”, 
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level and the 4.25-Mev level were enhanced some 
twenty times over the single-particle estimate. 
The spins, parities, and lifetimes of the first 
three levels of Ne*° suggest that they form the 
first three members of a perturbed rotational 
band. Preliminary evidence for a level in Ne”° 

at about 7.6 Mev has been obtained recently in a 
study of gamma-gamma coincidences" and a mag- 
netic spectrometer analysis’ of the alpha parti- 
cles from the reaction C!*(C™*, ~)Ne”°. This level 
has some of the properties required for a 6+ level 
but is shown dashed in Fig. 1 because of the pre- 
liminary nature of the evidence. 

Measurements?* have been made which showed 
that the spins and parities of the 5.63-Mev level 
and of a level reported recently by Adams et al.” 
at 5.80 Mev are 3- and1-, respectively. These 
assignments together with the measured values’® 
of f,/P of 0.07(+0.01) and 0.003 (0.003) show 
that both these levels, like the 2- level at 4.97 
Mev, have strongly inhibited E1 de-excitation 
widths. Bearing in mind the rotation-like band 
based on the ground state of Ne”°, it was conjec- 
tured that the 2-, 4.97-Mev level, the 3-, 5.63- 
Mev level, and the (5-), 8.84-Mev level’® formed 
a negative-parity band. A 4- level would then 
be expected at an excitation of about 7 Mev. The 
enhancement of £2 transitions within the ground- 
state rotational band and the strong inhibition of 
the £1 transitions from the 4.97-, the 5.63-, and 
5.80-Mev levels’®»*® suggested that the 4- level, 
if it exists and if it is indeed a member of the 
postulated rotational band, should decay by E2 
transitions to the 2-, 4.97-Mev and the 3-, 5.63- 
Mev levels in strong competition with a possible 
E1 decay to the 4+, 4.25-Mev level. 

A new level in Ne*° at 7.02 Mev, exhibiting sev- 
eral of the expected properties, has now been lo- 
cated using the C!#(C'*, wy)Ne”° reaction. The new 
level shows gamma-ray transitions to the 2-, 
4.97-Mev level and the 3-, 5.63-Mev level with 
relative intensity approximately 2:1. It is inter- 
esting to note that if these gamma-ray transitions 
are pure electric quadrupole the collective model 
predicts’” a relative intensity of 3:1. No other 
cascade gamma rays with intensities greater than 
15% of that for the transition to the 4.97-Mev lev- 
el have been observed. Alpha-gamma angular cor- 
relations are consistent with spin and parity of 
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4- but do not conclusively prove this assignment. 
Recent measurements” of alpha-particle spectra 
using a magnetic spectrometer at 0°, 10°, and 50° 
to the beam show that the 7.02-Mev state very 
probably has parity (-1)7 + 1 

In Fig. 2 the energy levels of Ne” are plotted 
as a function of J(J+1). In this plot one can see 
deviations from the simple rotational expression, 
E(J)=(#/23)J(J+1). The lines connecting states 
of each rotational band are expected to be approxi- 
mately parallel and this appears to be the case. 
The average slope of the lines, which is equal to 
i?/29, is approximately 150 kev. This may be 
compared with a slope of about 135 kev expected 
for the rotations of a rigid spheroid of distortion 
parameter’ 6~0.3. The dashed lines shown in Fig. 
2 are possible extensions of the rotational bands. 

The absence of a 1- level below the 2-, 4.97- 
Mev level indicates that the band based on that 
level can be assigned a K -quantum number of 2. 
Possible levels of spin and parity 2- and 4-, which 
might be members of the band based on the 1-, 
5.80-Mev level, have not yet been observed and 
the present evidence suggests that they do not 
exist. This point is of interest because the 1-, 
5.80-Mev level and the 3-, 7.19-Mev levels may 
then be the first two members of a band with K =0. 
If negative-parity bands with K =0 and K =2 can be 
established conclusively, then it is possible that 
these bands may be explained as rotations of the 


Ne”® nucleus which is undergoing surface octu- 
pole vibrations.’® The alternative explanation 
that the negative-parity bands are based upon 
states of particle excitation must also be con- 
sidered. According to the Nilsson mbodel,?* the 
lowest of these is expected to be a configuration 
in which one particle is removed from the [101] 
level and put into the [211] level. Such particle 
configurations would give rise to two bands of 
K=1-, J=1, 2, 3, etc., andK=2-, J=2, 3, 4, 
etc. Calculations” indicate that these bands 
would not occur as low as 5 Mev in Ne” and the 
absence of the J=2- level at ~6.3 Mev argues 
against a K=1 assignment to the band based on 
the J=1- level at 5.8 Mev. 

Clearly more experimental evidence is re- 
quired on the existence of other band members 
and their properties before a complete account 
of Ne”° is possible. However, the collective 
model gives a good account of the properties of 
the low-lying levels of Ne*° which are known at 
present. 

We would like to acknowledge the stimulating 
interest shown in this work by Dr. L. G. Elliott. 
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PARTICLE PRODUCTION BY 10-30 Bev PROTONS INCIDENT ON Al AND Bet 


W. F. Baker, R. L. Cool, E. W. Jenkins, T. F. Kycia, S. J. Lindenbaum, W. A. Love, D. Liiers,* 
J. A. Niederer, S. Ozaki, A. L. Read, J. J. Russell, and L. C. L. Yuan 
Brookhaven National Laboratory, Upton, New York 
(Received June 21, 1961; revised manuscript received July 12, 1961) 


High-energy secondary particles produced by 
the nuclear interactions of 10-30 Bev protons 
have been investigated at the Brookhaven alter- 
nating gradient synchrotron (AGS).! This experi- 
ment was carried out with aluminum and berylli- 
um internal targets (according to availability) at 
production angles? of 43°, 9°, 13°, and 20°. The 
momentum of the particles in a beam was selected 
by magnetic deflection and the velocities deter- 
mined with focusing, gas Cerenkov counters.® 
The momentum acceptance and angular divergence 
of the beams were fixed by scintillation counter 
telescopes. Typical values of these parameters 
were Ap/p~+2% and AQ~2x10~ steradian. The 
ratio of quadruple coincidences (S,S,S,C) to triple 
coincidences (S,S,S,) corrected when necessary 
for Cerenkov counter inefficiency gave the rela- 
tive abundance of the particle of selected mass in 
a beam. 
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Most of the data were taken with the machine 
operated at 29.5+ 0.5 Bev with an average intensi- 
ty of ~4x10"° protons per pulse accelerated each 
2.4 seconds. Measurements were also made at 
AGS energies of 25, 20, and 10 Bev, at which en- 
ergies the machine can be operated at 2.0-, 1.6-, 
and 0.8-second repetition periods, respectively. 

Figure 1 shows the momentum distributions of 
protons, 7* and 1~ mesons, corrected for 7 decay 
between the target and the final scintillator in or- 
der to represent the pion-production rates at the 
target. A rough correction has been made for the 
number of decay muons counted as pions. This 
correction was less than 5% of the measured pion 
intensity for beams of momentum >4 Bev/c. The 
muon contamination increased rapidly with de- 
creasing momentum below 4 Bev/c; at 2.5 Bev/c 
the contamination was ~12% and at 1.25 Bev/c it 
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FIG, 1. The momentum spectra, (a) of 7* mesons and protons for 30-Bev incident proton energy, and (b) of 7* 
mesons for 10- and 20-Bev incident proton energies. The ordinate is the number of particles per steradian per 
Bev/c per circulating proton. The dashed portions of the 4#° spectra indicate the regions where the correction 
due to fringing field effects exceeds 15% of the measured value of d*N/dpdQ. 
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ber of circulating protons in the AGS which has 
been calibrated by electrical methods to an ac- 


mately represents the number of particles pro- 
duced per inelastic interaction with the target 
nucleus per steradian per Bev/c. Fluctuations 
in target efficiency may have been as much as 
+10%. 


of the AGS strongly influences only the trajectory 
of the nominal 43° beam. Thus in this beam the 
actual production angle varies from 33° for nega- 
tive particles of momentum 6 Bev/c to 53° for 


duction angle is 4.7° for very high momentum (i.e 
~20 Bev/c) particles of either sign of charge.” 
The fringing field of the AGS and the beam colli- 
mator combine to reduce the effective solid angle 
seen by the counter telescope. Corrections for 
these effects have been applied to the measured 
data; the dashed portions of the 43° curves indi- 













each target of ~50%. Thus, 2xd*N/dpdQ approxi- 


Of the four beams measured, the magnetic field 


positive particles of momentum 2 Bev/c; the pro- 






was ~40%. The curves are normalized to the num- 
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cate the regions where these corrections exceeded 
15%. Figures 2 and 3 and Table I show the ratios 
of the intensities Kt/7*, K~/n~, and p/n~, meas- 
ured in each beam at various momenta. The ratios 
measured experimentally have been corrected for 
K and 7 decay between the target and the final scin- 
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10 


FIG, 2. The production ratios at the target of 
K*/n*, K~/x~, and p/1~ at various momenta, at 43° 
production angle. The ratios have been corrected 
for K and 7 decay between the target and the final 
scintillator and for muon contamination. 


Table I. The production ratios at the target of K*/ a*, K~/x~, and p/n™ at various momenta for production angles 
of 13° and 20°. The AGS internal energy =30 Bev. The ratios have been corrected for K and 7 decay between the 
target and the final scintillator and for muon contamination. 
























Production Momentum Par Production ratio, at target 
angle Target (Bev/c) K /n K /n p/n 

13° Al 3.0 0.27 +0. 02 
4.5 0.013 +0.002 
6.0 0.29 +0.05 0.062 +0.010 0.014 +0. 002 
9.0 0.31 +0. 04 0.048 +0.010 0.013 +0. 007 

20° Be 1.25 0.027 +0.012 0.012 +0. 005 
2.5 0.230 +0.012 0.088 +0. 009 0.0051 +0. 0015 
3.0 0.088 +0.010 
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FIG. 3. The production ratios at the target of K*/r*, 
K~/x~, and p/x at various momenta, at 9° production 
angle. The ratios have been corrected for K and 7 
decay between the target and the final scintillator and 
for muon contamination. 


tillator so that the values plotted represent the 
production ratios at the target. 

It can be seen from Fig. 1 that, where compari- 
sons can be made, the characteristics of the pion 
spectra from the beryllium target and from the 
aluminum target are almost the same, except 


that 
( dN ) ik eu) 
a= oma ~~ 1.1 -1.4., 
dpan) ,, | \dpdo) ,, 


Part of these apparent differences between Al and 
Be targets may be due to relative monitoring and 
beam normalization errors and to fluctuations in 
target efficiency. The ratios K*/n*, K~/a~, and 
p/7~ can be seen from Figs. 2 and 3 to be insensi- 
tive to the difference between an aluminum and a 
beryllium target. In view of this insensitivity of 
particle production to target material and from 
general considerations based on the optical model, 
it appears that an analysis of the data in terms of 


a collision of the incident proton with a single nu- 
cleon of the target nucleus may be of interest and 
may reveal some significant characteristics of 
high-energy proton-nucleon interactions. 

If one transforms the pion spectra to the proton- 
nucleon c.m. system, one finds that the pion spec- 
tra exhibit a low-energy peak $500 Mev/c in the 
c.m. system, although the available kinetic energy 
<6 Bev. Furthermore the angular distribution is 
peaked in the forward direction. The c.m. angular 
and momentum distributions of the pions appear to 
be consistent with a limited mean transverse mo- 
mentum transfer £0.5-1 Bev/c. 

From Figs. 2 and 3 it is clear that the K*/z* 
ratios are not strongly sensitive to incident pro- 
ton energy and tend to increase toward higher Kt 
momenta. For 20-30 Bev incident proton energy, 
typical measured values of the K*/n* ratio vary 
from 0.05 to 0.3. On the other hand, the K~/n~ 
ratio falls rapidly with decreasing incident proton 
energy and also falls rapidly with increasing K~ 
momentum. Typical measured values of the 
K~/n™~ ratio vary from 0.08 to 0.02. At a given 
production angle the K~/n™ and K*/a* ratios tend 
to approach equality at very low momenta. These 
results may imply that the production of low-mo- 
mentum K mesons occurs predominantly through 
K-pair production, while in the production of high- 
momentum K* mesons associated production 
(K -hyperon production) predominates. 

For p/1~ ratios typical values are ~0.01-0.015 
at low p momenta, falling to ~0.001 at high mo- 
menta. The p/m™~ ratios are sensitive to the in- 
cident proton energy, dropping rapidly at the 
higher ) momenta when the incident proton energy 
is reduced from 30 Bev to 20 Bev. This result is 
understandable when one recalls that the p-p 
threshold in the c.m. system is 1.88 Bev. 

Comparisons can be made between the results 
presented here and the results of some investi- 
gations’ performed at CERN at 25-Bev proton en- 
ergy using an Al target. At 6° lab production angle, 
three CERN K~/z” ratios were reported —at 5, 6, 
and 8 Bev/c. These data seem to lie reasonably 
between our data at 30-Bev and 20-Bev incident 
proton energy. Of the four CERN K*/n* ratios 
reported by von Dardel et al.,’ the ratio at 5 Bev/c 
is higher than our result by a factor ~1.5. How- 
ever, the CERN K*/n* ratios at 6 and 8 Bev/c 
agree reasonably well with our results. 

The CERN p/z~ ratios at 6° lab production angle, 
at 6, 8, and 11 Bev/c, are in agreement with our 
data, but the CERN point at 16 Bev/c appears to 
be higher than the trend of our data. The p/7~ 
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ratios measured at CERN at 16° lab production 
angle by Cocconi et al.’ are reasonably consistent 
with the trend of our measurements at production 
angles of 13° and 20°, 

A more complete account of this work will be 
published elsewhere. 

We wish to acknowledge the invaluable coopera- 
tion of the AGS staff and operating crew which 
made this investigation possible. 
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FAST ATOMIC TRANSITIONS WITHIN » -MESONIC HYPERFINE DOUBLETS, 
AND OBSERVABLE EFFECTS OF THE SPIN DEPENDENCE OF MUON ABSORPTION* 


R. Winston and V. L. Telegdi 
The Enrico Fermi Institute for Nuclear Studies, The University of Chicago, Chicago, Illinois 
(Received June 19, 1961) 


It has been pointed out by Bernstein, Lee, Yang, 
and Primakoff,’ that the two members F,, =/+ of 
the hyperfine doublet ground state of a u-mesonic 
atom with nonzero nuclear spin J could exhibit 
different disappearance rates A + aS a consequence 
of the possible spin dependence of the interaction 
responsible for muon absorption. In particular, 
the now currently favored “universal” V-A inter - 
action would lead to A,<A_. Bernstein et al.’ 
made the implicit assumption that atomic proc- 
esses inducing transitions between the two states 
F, were negligibly slow. With this assumption, 
they predicted that a logarithmic plot of the time 
dependence of the electron rates from negative 
muons stopped in a monoisotopic target of non- 
zero J should exhibit a positive curvature, cor- 
responding to the superposition of the two expo- 
nentials characterized by A, and A_. Experi- 
mental detection of such a curvature could serve 
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as a proof of the spin dependence of the absorp- 
tion interactions but would give no clue to its 
more detailed nature. 

It was subsequently pointed out? that, at least 
for the targets of experimental interest, the as- 
sumption of a negligibly slow transition rate be- 
tween the two states F, is not a valid one. It was 
shown that the magnetic interaction between the 
“core” (=nucleus + muon) of the mesonic atom 
and its outer s electrons provides a sizeable con- 
version rate R between those two states. An ex- 
plicit calculation for the case of ,,Al”’ was pre- 
sented in reference 2, made under the assumption 
that only the 3s (conduction) electrons of this tar- 
get would effectively contribute to R. The magni- 
tude of R so calculated was comparable to that of 
the difference |A_-A,|=lA4A\l. As was empha- 
sized in reference 2, this parent-daughter rela- 
tionship would lead to a negative curvature in the 
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experiment proposed by Bernstein et al. if the 
state F_ disappeared faster, which is precisely 
the prediction of the V-A interaction for a nucleus 
with positive magnetic moment yu r 

Thus, the curvature experiment could possibly 
provide the as yet missing evidence that the V-A 
interaction governs muon capture. Some prelim- 
inary experimental evidence from Chicago, indi- 
cating the actual observation of negative curva- 
ture in Al and P, was discussed at the 1960 Roch- 
ester Conference.* Subsequent more refined meas- 
urements* performed in this laboratory indicate, 
however, that this preliminary “evidence” was 
due to spurious instrumental effects, the nature of 
which is now understood. This situation leads 
to the following questions: (a) Granting that muon 
capture is governed by the V-A interaction, what 
modification of the plausible rate calculation in 
reference 2 can account for the nonobservation of 
negative curvature? (b) Granting that a modified 
(and presumably reliable) rate calculation can ac- 
count for this, what modification of the experi- 
ment originally proposed by reference 1 would 
lead to an observable negative curvature, thus 
confirming V-A? The present note purports to 
answer these questions. 

(a) Conversion rate calculation. The problem 
of the transition rate between the states F,, of in- 
terest, separated by an energy €, is entirely anal- 
ogous to that of an M1 transition between two nu- 
clear states, the “core” playing here the role of 
the nucleus. However, the radiative transition 
rate R’ is here readily calculable’®: 


R’'= ffm I/(2I + 1)e°, (1) 


while the conversion coefficients (that would yield 
R) are not available for the very soft («<1 kev) 
photons in question from standard tables.®° We 
must, as was done in reference 2, estimate the 
rate R; . ¢ at which the magnetic hyperfine inter- 
action between the “core” and an electron of the 
mesonic atom leads to electron ejection: 


=4ra“R'e]u (0) 1" lu, (0) 1?o,(k), 
(2a) 





= >ks 


R 
NDyj_g> kp 
= “lpi. > = - 2 
ira Re °(33/2)I J, May” Upp” ar Pp (R), 
(2b) 
R 
"Py. > kp 


=ly,-lp:-73 ar -3 2 
sta *R’€~°(69/10)1 f, Uap” upp” ari Ppl). 
(2c) 





The approximation in the above expressions con- 
sists in taking the proper nonrelativistic limits, 
neglecting retardation effects and dropping con- 
tributions from nuclear magnetic moments. The 
rates in (2) are already summed over the elec- 
trons of the pertinent atomic shells, assumed to 
be filled; the total conversion rate is clearly 


R=DR (3) 


nlj +klj” 


where the sum is taken over all ionizable shells 
(nlj); ionizable means that the “edge” E, ,. lies 
below the available energy 


€=E(F ,) -E(F _) 
~*m yp 7/4 ptt DMZ es (4) 


nlj 


~ 2,2 
$a *, 
where Zor’ =Zogs*/Z taken from yp capture.” The 
replacement of Z* by Zg¢,* corrects for the finite 
nuclear size. For a completely general discus- 
sion, see Bohr and Weisskopf.® 
We estimate R,, the conversion rate for a given 
shell n, by the following approximate expressions 
valid for |€-V,|«Z"/2: 


lu ,(0)1?*42"(Z" -S)?/n°, Ia, (0) 7p ,(k) 42", 


(5a) 


where Z’=Z -1, the charge of the mesonic atom 
core, and 


°o 


=3 2 ~7'2(7' _o\2 3 
If Uys Upptarl Pr (k)=Z'*(Z' -S)?/9n*, (5b) 
so that 
z., es » ,,, aiead! (5c) 


Here, S=inner, V,=outer Slater screening con- 
stants. These estimates are straightforward gen- 
eralizations to off-diagonal matrix elements of 
Goudsmit’s semiclassical hyperfine formula, jus- 
tified by Fermi and Segré via the JWKB method.® 
For the shielding of the emitted s electron, the 
remarks of Rose’ are pertinent. In reference 2, 
R was calculated from (2a), assuming that only 
the 3s electrons of Al (magnetic My shell conver- 
sion) could be ejected. Inasmuch as €(Al) ~260 ev 
while E255 = 88.60 ev (exp. L; - edge of Mg”), 
this assumption was in error. Using (3), (4), and 
(5), one finds: 


R,~3.2 x10" sec™ for ,,Al” (6) 
(Z'=12, S=4.15, V,=120 ev). 


The numerical value (6) far exceeds that calcula- 
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ted in reference 2, viz., R=6.2x10° sec™. It is 
two orders of magnitude larger’ than |AA| and 
this fact, as we shall show later, accounts for the 
absence of observed negative curvature in Al. For 
sP™ the situation is entirely similar. 

R> |AA\ is not restricted to Al and P but holds 
quite generally for odd-Z nuclei of interest through- 
out the periodic table. This is due to the fact that 
whenever s conversion can occur, it does so ata 
rate that greatly exceeds |AA| (specifically dis- 
regarding nuclei with anomalously small yu p- 

Since R,, increases with Z more rapidly than does 
IAAI ~Zage’, it is sufficient to show this for the 
lowest Z at which a particular n,s shell can first 
be ionized. Thus, for Z<9, conversion occurs 
first at Z +5; for ,B“, R,~1.7x10* sec™, AA 
~5x10° sec. For 9<Z<15, 2s conversion al- 
ways occurs; for ,F'®, Ry, ~2.9x10° sec”, AA 
~1.1x10° sec™. For 17<Z<39, 3s conversion 
always occurs; for ,,CI®, Rg,~7.5x10° sec™, 
AA~2.3x10° sec™. For 39<Z<75, 4s conver- 
sion always occurs; for ,,Y°®, R,~4.110" sec™, 
AA™~2.8x10° sec™. Since the observable effects 
will be shown to decrease as 1/Z, we do not con- 
sider very high values of Z. 

The high conversion rates (R > |AA|) discussed 
exclude practically observable effects in the time 
distribution of decay events, quite independently 
of the nature of the y-capture interaction. Assum- 
ing uy>0, the time dependence of the populations 
in the F, states is governed by the equations 


asihenmpe oO. 


n_(t) =[n_(0) +n,(0)R(R -AA) Je 3! 
-1 -(A,+R)t 
-n,(O)R(R-AA) e mee » (7) 
where n,(0) are the statistical weights of the states 
F. In the limit that R> IAAI, the time distribu- 
tion of the decay events (i.e., the electron rate) 
to order AA/R is 
d -At - 
No“ (t) ~e A. (1 -Ae ~ (8) 
where A=n,(0)AA/R. The smallness of the ampli- 
tude A together with its rapid decay ( in ,,Al*’, A 
~2x107°, R“*~3 x10™ sec) accounts for the ab- 
sence of observed negative curvature in ,,Al”" and 
isP". As emphasized above, this smallness is 
not restricted to these nuclides. 
(b) Experiments leading to observable negative 
curvature. Contrary to the study of the time dis- 
tribution of the decay events, as originally pro- 





106 


posed in reference 1, a study of the time distri- 
bution of capture events (capture neutrons and 
gamma rays) can lead to practically observable 
spin-dependence effects and, in particular, ex- 
hibit negative curvature. One has for these events 
from (7), again in the limit R > AA, 


Rt 


nig ~e A“ -a *, (9) 


where A’=[n,(0)AA/A@P]/1 +n, (0)AA/A“P}, 
ACP en ,(0)A,°?P +n_(0)A_©9P = mean capture 
rate. The striking difference between (8) and (9) 
is that A’ may be a large amplitude of order unity, 
It is easy to see how this comes about: 

When |AA|/R<1, the muon population is for 
t>>R~ captured almost entirely from the F_ 
state at a rate A_©?Pe™ Extrapolating back 
to t=0, this would require an initial capture rate 
a. whereas the true initial capture rate is 
ACP, Therefore, 1-A’=A°@P/A_©4P, where A’ 
=[n,(0)AA/AC@P]/[1 +n.(0)AA/AC@P] as in (9). 
The same argument applied to decay events shows 
that A = 0 to zeroth order in AA/R, since the de- 
cay probability is the same in both F states. 
Since A’ decreases as 1/Z, the effect is best 
studied in low-Z elements. For statistical ef- 
ficiency, R should not greatly exceed A_. Several 
elements (e.g., ,F*®) appear to be suited in prac- 
tice for such experiments. 

We conclude with two remarks: 

(i) If one had A,/R >1, targets with /=3 (e.g., 
F*® and ,,P™) should exhibit** a u-decay asym- 
metry parameter a(J) half as large as a(0), the 
corresponding parameter for J/=0. With A ,/R =1, 
one should see a “damping” of the precession in 
the F, state, while A,/R <1 would imply a(3)=0. 
In the case of ,,P** (A,/R~0.03), both a(3)/a(0) 
=0.5** and a(4)/a(0) = 02° have been reported, 
while recent Liverpool results on ,F®(A,/R = 0.25) 
are consistent with a(4) /a(0) = 0.7* 

(ii) With A,/R «1, virtually all muons should 
be absorbed from the F _ state (assuming yu r 0). 
The V-A interaction predicts A_>A,; thus ona 
plot of the reduced capture rates A°@P/Z a¢¢* vs 
neutron excess (A-Z)/2A, such nuclei as ,F’’, 
Na’, ,,Al’’, and ,,P™ should lie above the “Prim 
akoff line.”“” Some evidence for this has been ob- 
tained recently,* but a quantitative comparison 
with theory does not appear possible at the prese 
time. 

We wish to thank Professor H. Primakoff and 
Professor G. Wentzel as well as J. L. Lathrop 
and R. A. Lundy for helpful comments and dis- 
cussions. 
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MEASUREMENTS OF MUON DISAPPEARANCE RATES vs TIME IN C, Mg, Al, Si, AND PT 
J. L. Lathrop, R. A. Lundy, V. L. Telegdi, R. Winston, and D. D. Yovanovitch 


The Enrico Fermi Institute for Nuclear Studies, The University of Chicago, Chicago, Illinois 
(Received July 13, 1961) 


At the 1960 Rochester Conference we reported! 
some preliminary evidence for the observation of 
a hyperfine structure (i.e., spin dependence) ef- 
fect in the capture of muons by ,,Al”” and ,,P™*. 
This effect consists in a (negative) curvature® 
K of the plot of f(t), the logarithm of the electron 
rate from these targets, vs time, 


K =K(0) =f"(1+f) ™ at t=0, (1) 
and depends on the physical rates as follows: 
K=K(AA, R) = -27*"n,(AA/A)(R-n_AA)/A, (2) 


where m, are the statistical weights of the hyper- 
fine states F, =J+3, A=n,A,+n_A_ the mean dis- 
appearance rate, AA=A_-A,, and R is the con- 
version rate from F, to F_ (assuming a positive 
nuclear magnetic moment).* 

We have now repeated these “curvature” meas- 
urements under greatly improved conditions, and 
have found that, to within experimental error, K 
=0 for Al and Mg. We used the experimental set- 
up schematically indicated in Fig. 1. The most 
essential improvements were the following ones: 

(a) Comparison of the Al target (J= §) data with 


(e-)=5641 
(u-)=2345 


tv 


Mg data (mostly J/=0) obtained under essentially 
identical conditions. In the past, C reference 





Inches 


FIG. 1. Experimental setup. 1, 2, 3, 4,5, and6, 


square scintillation counters; A; = j-in. Cu moderator, 
A, =2-in. graphite moderator, A;=4-in. polythene ab- 
sorber. Pb=Pb collimators. (1456) =electron coinci- 
dence, (2345) = muon coincidence. 
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Table I. Experimental results on curvature in Mg and Al. 
Accidental background Carbon background 
Target Kx10° (% at t=0) (% at t =0) 
Mg +(1,1+1.0) 0.209 +0.002 0.09 +0.38 
Al +(2.4 41.7)" 0.295 +0.004 0.03 +0.29 





*To be compared with the value K=-(4.4 +1.0) x1073 given in reference 1. 


targets had been used to get “uncurved” data; 
however, C targets can show no background from 
muon decays in carbon-—an essentially inevitable 
phenomenon with all other targets. This carbon 
background could simulate negative curvature 
through reduced reliability in fitting the time 
distribution. 

(b) Improved time range of observation. To get 
a better evaluation of the flat, accidental back- 
ground, an interval of 9 usec was observed at 
t<0 (i.e., before the arrival of each muon). A 
longer (16-sec) time range at t>0 was made 
available for observation in order to make the 
data statistically more valuable for separation 
(by least-squares program) into flat accidental 
background and exponential carbon background. 
Accidental backgrounds from ¢<0 data and from 
fit to t>0 data were required to match. 

(c) Improvements in “digitron,”*** the clock used 
to measure the distribution of uw -e time intervals. 
These improvements, to be described in detail 
elsewhere,° consisted primarily in replacing the 
pulsed oscillator with a continuous oscillator 
and scaler controlled gate combination, and in 
general using dc coupled transistorized circuits. 

(d) Reduction of accidental background, achieved 
both by providing additional protective circuits in 
“digitron” and by using a high duty-cycle (~20%) 
muon beam.® 

(e) Reduction and direct measurement of carbon 
background. Directly measured carbon background 
and carbon background from least-squares fits to 
t>0O data were required to agree. Reduction was 
acheived by covering scintillators 4 and 5 with thin 
(1-mil) Ag foils and not using Scotch tape in addi- 
tion thereto, as well as by increasing their anti- 
coincidence efficiencies. The short-lived muons 
stopping in these Ag foils essentially do not con- 
tribute electrons to the ¢>0 range used in data re- 
duction. The direct carbon background measure- 
ments were made with an “equivalent” Cu target. 
The short (Cu) and the long (C) components sepa- 
rate easily. Two layers of Scotch tape (16 mg/ 
cm’), added on the upstream face of counter 5, 
could readily be detected in the presence of a 9 
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g/cm? Cu target. This purposely created carbon 
background could be used to gauge the carbon 
background present in the actual (Mg and Al) r 
To make the “curvature” (Al) and the “control” 
(Mg) runs as closely comparable as possible, the 
Al target consisted of 5 equally spaced j-in. thi 
sheets filling the same volume (6 in. x6 in. x2 in. 
as the Mg target. Table I summarizes our most 
reliable data obtained with these targets; to withi 
experimental error, K =0 in both Al and sign | 
contrast to our preliminary evidence for Al. A 
recent calculation of the conversion rate’ yields 
R(Al) =3.2x10" sec™. With this large value of R 
K(t) becomes a rapidly varying (decreasing) func 
tion of time.* Under these circumstances the 
“curvature” K obtained by fitting the electron 
time distribution over several mean lives is no 
longer given by (2) but rather by K = -2™*n, 
x(AA/R)(A/R). With the above value of R and 
the estimate? AA/A=0.17, one has K~ -5x107. 
This is in agreement with the null result in Tabl 
I. 










In the course of our search for spin-dependent 
effects in muon capture, we have obtained rather 
precise pu lifetimes for C, Mg, Al, Si, and P. 
These, as well as some quantities derived from 
them, are summarized in Table II. 


Table II, Precise u™ lifetimes for C, Mg, AI, Si, 
and P, 








Capture rate” 

Lifetime AAP x1 975 P 

Target (usec) (sec™*) A°*P ip) 
Cc 2.041 +0.005 0.361+0.01 17.8+0.5 
Mg 1.071 +0.002 4.80 +0.02 19.4+0.1 
Al 0.864 +0.002 6.62 +0.03 21.4+0.1 
Si 0.767 +0.002 8.50 +0.03 19.9+0.1 
Pp 0.635 +0.002 11.21 +0.05 21.0+0.1 





* Statistical error computed from least-squares fit 
to, data. 
Computed using the free muon decay rate, Agec 
= (4.539 0.005) x10® sec, 
(p) as defined and given by J. C. Sens, Phys. Re 
113, 679 (1959). 
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SEARCH FOR HIGH-ENERGY COSMIC GAMMA RAYS* 


Thomas L. Clinet 
Department of Physics and Laboratory for Nuclear Science, 
Massachusetts Institute of Technology, Cambridge, Massachusetts 
(Received June 19, 1961) 


This note describes an attempt to detect gamma 
rays of cosmic origin in the energy region appro- 
priate to 7° decay. Since 7° mesons are produced 
in nucleon-antinucleon annihilation, the existence 
or nonexistence of a gamma-ray flux from certain 
portions of the sky bears upon questions such as 
possible collisions between galaxies and antigal- 
axies? and cosmological models which postulate 
matter and antimatter creation.” New upper lim- 
its are set on the creation rate and on the density 
of interstellar antinucleons. 

The existence of high-energy gamma rays in 
the primary cosmic radiation was first investi- 
gated by Schein, Jesse, and Wollan using G-M 
tubes,® by Hulsizer and Rossi using ionization 
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FIG, 1. Cross section of the 
apparatus. The geometrical 
figure for this telescope was 
1.7 cm? sr. 
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chambers,‘ and by Critchfield, Ney, and Oleksa 
using cloud chambers.® These experiments set 
an upper limit for the flux of the electron-photon 
component above about 1 Bev at about 1% of the 
primary cosmic-ray flux. G-M telescopes were 
used by Perlow and Kissinger,® and more recently 
by Danielson.” All of these experiments suffered 
the disadvantages of either a high proportion of 
locally-produced background or an energy sensi- 
tivity which did not extend significantly into the 
70-Mev region, and none had a directional sur- 
vey of the sky as its purpose. 

A cross section of the balloon-borne detector 
that was used in the present experiment is shown 
in Fig. 1. Incoming gamma rays, collimated by 
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the lead shield, produced electron pairs in either 
of two removable mercury radiators. The upper 
radiator, when inserted in flight, was used to 
measure the attenuation of the beam in order to 
identify the measured flux as photons. The lower 
radiator produced the pairs that were detected 
by the scintillator -Cerenkov combination at the 
bottom. A pulse corresponding to the traversal 
of two minimum -ionizing particles in the cesium 
iodide crystal was required in coincidence with 
a pulse from the Lucite Cerenkov counter. The 
pair detector was surrounded by an anticoinci- 
dence scintillator inside the lead shield. This 
anticoincidence requirement excluded incident 
charged particles as well as electron pairs or 
showers projected upwards from the detector. 

In addition, a large anticoincidence scintillator 
outside the lead shield was used to reduce the 
background produced by cosmic-ray interactions 
in that shield. The apparatus was rotated in 
flight to various zenith angles in order to pour 
the mercury in and out, to measure the zenith- 
angle dependence of the gamma-ray intensity 
and to measure the gamma-ray albedo. 

The detection efficiency of the instrument was 
estimated as a function of gamma-ray energy 
from the known pair-production cross sections, 
electron multiple-scattering expressions, etc. 
This curve was combined with the experimentally 
determined energy spectrum of pair-producing 
gamma rays near the top of the atmosphere as 


measured in the emulsion experiments of Carlson, 


Hooper, and King,’ and of Svensson.® The result 
was an effective efficiency of about 0.17. This 
efficiency differed by less than 10% with the sub- 
stitution of the gamma-ray spectrum from low- 
energy proton-antiproton annihilation. The ade- 
quacy of the efficiency estimate was confirmed 
and the angular response of the instrument was 
measured in a 7° decay gamma-ray beam at the 
MIT synchrotron. 

The detector was balloon-borne at 55° geomag- 
netic latitude on July 1, 1960, for about one day of 
which 12 hours were spent at 8.5 g cm~™ atmos- 
pheric depth. Evidence that gamma rays were 
detected is as follows: The counting rate as a 
function of atmospheric depth, shown in Fig. 2, 
followed a transition curve having a maximum 
near 180 g cm™, entirely similar to that of the 
electron-photon component. The ratio of the 
counting rate with the upper mercury absorber 
inserted to the rate with this absorber removed 
was 0.40+0.07, while the expected ratio for pure 
gamma-ray detection was 0.35. The gamma-ray 
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albedo flux was measured and found to be in agree- 
ment with the results of the emulsion experiment 
of Svensson, having a value of about half that of 
the vertically downward gamma-ray flux at 14 g 
cm™, 

The gamma-ray intensity versus atmospheric 
depth in the region 0 to 40 g cm™ is shown in the 
inset of Fig. 2. The four measurements with the 
smallest indicated errors were taken at four ze- 
nith angles while the balloon floated at 8.5 g cm™, 
Since the corresponding atmospheric depths are 
small compared with any of the characteristic 
interaction lengths involved, a linear extrapola- 
tion to 0 g cm™ is justified. The extrapolated 
intensity, determined by a least-squares fit, is 
(1+3)x107> cm™ sec sr™ and the correspond- 
ing upper limit to the gamma-ray intensity inci- 
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FIG. 2. Gamma-ray intensity versus atmospheric 
depth. (The axis of the detector was clamped at 30- 
degree zenith angle throughout the ascent, hence the 
15% correction to depth.) The insert shows a linear 
plot of the intensity for small atmospheric depths; the 
units are the same for both graphs. 
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dent upon the top of the atmosphere, for a 95% 
statistical confidence limit, is 7x107° cm~* sec™ 
sr. This value must be regarded as an average 
over the directions scanned, i.e., as a limit to 
the essentially isotropic intensity from the north- 
ern celestial hemisphere. 

From a sidereal arrival direction analysis, the 
incident flux from the direction of Cygnus-A, hav- 
ing subtracted off the secondary flux, was found 
to be (-2+7)x10™“ cm™ sec™, with the corre- 
sponding upper limit of 1.2 x 107° cm~ sec™ for 
a 95% statistical confidence limit. A previous 
upper limit, based upon the detection of one gam- 
ma ray, was set at 5x10~° cm™ sec™ by Brac- 
cesi, Ceccarelli, and Salandin,’® using emulsions. 
This corresponds to an upper limit of 1.5107 
em™~™ sec™ for a 95% statistical confidence limit. 
Estimates of the flux from Cygnus-A, based on a 
galaxy-antigalaxy collision model, range from 
Morrison’s value of about 1 cm™ sec™ to Save- 
doff’s value," which is yet three orders of mag- 
nitude below the upper limit set by this experi- 
ment. 

Proton-antiproton annihilation taking place any- 
where in the cone defined by the aperture of the 
instrument could also contribute to a gamma-ray 
intensity. The lifetime of an antiproton against 
annihilation with interstellar protons can be taken 
as T=[nv(o,c/v)]“, where n is the proton number 
density, v is the relative velocity, and o=0,(c/v) 
is the annihilation cross section, in which o, is 
about 10-8 cm?. This lifetime is much greater 
than the time scale of the universe (1.4x10"° 
years) for proton densities appropriate to inter- 
galactic space (n<10~* cm™), and is much less 
than the time scale of the universe for proton 


densities appropriate to the galaxy (7 > 107? cm“). 


Within wide limits of proton densities, the annihi- 
lation frequency would therefore be constant with- 
in the galactic volume, so that the resulting anni- 
hilation gamma-ray flux depends only upon the 
geometrical shape of the galactic portion sub- 
tended by the aperture of the detector. The anni- 
hilation frequency, S, and the gamma-ray inten- 
sity, J, are related by S=4nJ(mRf)™, where m 
23 is the average number of gamma rays per 
annihilation, R is the distance to the surface of 
the galactic disk in the direction scanned, and f 

is a geometrical factor which takes into account 
the shape of the galaxy within the cone scanned. 
For the direction perpendicular to the galactic 
Plane, R =5x10"° cm and f=1; for the direction 

of the galactic anticenter, R =2x10”* cm and 
f=0.15. 


During the time that the balloon was near max- 
imum altitude the aperture of the instrument 
scanned the galactic anticenter at certain known 
times, but, due to a partial loss of azimuth in- 
formation, it was not possible to separate all of 
the data from the direction perpendicular to the 
galactic plane. However, an upper limit for each 
intensity can be established: The greater inten- 
sity would come from the galactic anticenter, so 
that a conservatively set upper limit results if 
the value for the isotropic intensity is treated as 
a value for the direction perpendicular to the 
galactic plane. With J<7x107° cm™ sec™ sr“, 
this gives S<6x107** cm™ sec™ as the upper 
limit to the annihilation frequency for the portion 
of the galaxy in the direction of the North Galactic 
Pole. This limit is independent of knowledge of 
the shape of the galaxy other than the perpendicu- 
lar distance to the surface. An upper limit to the 
gamma-ray intensity from the direction of the 
galactic anticenter, for a 95% statistical confi- 
dence limit, was calculated to be 1.8107? cm~ 
sec! sr“. This value of J gives S <2.5x107** 
cm™~ sec™ as the upper limit to the annihilation 
frequency for the galactic anticenter. 

The product of the annihilation frequency, S, 
and the antiproton lifetime, 7, is the correspond- 
ing upper limit to the antiproton number density, 
n. The result, which is a function of the proton 
number density, can be expressed as 7 =S(o,cn)™ 
<(8x10~°)/m cm. This upper limit to the inter- 
stellar antiproton density is a factor of ten below 
previous limits that assume n=1 cm™~.? 

The upper limit for the annihilation frequency 
of antiprotons can also be taken to be an upper 
limit for the creation rate of antiprotons. As 
long as the antiproton lifetime within the galaxy 
is short compared with the age of the galaxy, 
which seems certain for any proton density in 
excess of 10 em™, an equilibrium between 
annihilation and any postulated long time-scale 
production mechanism should be established. 

The proton production rate has been required by 
steady-state cosmology to be about 3 x10-*? cm™ 
sec™,? which is a factor of ten above the upper 
limit of 2.510773 cm~ sec™ for antiproton anni- 
hilation (or production) reported here. The value 
of this required production rate in steady-state 
theory, however, depends upon the experimental 
determination of, e.g., the average universal 
mass density, which at present is uncertain and 
may be as much as a factor of 10° below the value 
of 10-" to 10-?° g cm™ assumed earlier by cos- 
mologists.’” 
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SELF-CONSISTENT CALCULATION OF THE MASS AND WIDTH OF THE J=1, T=1, 17 RESONANCE 


Fredrik Zachariasen* 


California Institute of Technology, Pasadena, California 
(Received June 28, 1961) 


The existence of a J=1, T=1 resonance in 17 
scattering now seems to be very probable, though 
there are still rather large uncertainties in its 
position and width.’ From a theoretical point of 
view, one may believe this resonance to be due 
to the existence of an unstable vector meson,’ 
called p, or one may say the resonance is dy- 
namical.* If the p meson has an infinite bare 
mass, there is essentially no difference in the 
results of these two points of view.**® The princi- 
pal distinction between the two philosophies comes 
from the fact that if the p meson is really a new 
particle, its mass and coupling constant must be 
thought of as new independent parameters which 
can be chosen to be anything, while if the p is a 
dynamical resonance, the mass and coupling con- 
stant are determined by other parameters already 
in the theory. 

Within the “dynamical” philosophy, attempts have 
been made to calculate the mass and coupling con- 
stant of the p meson.**® These have not been en- 
tirely successful; furthermore, the most complete 
and careful attempt to do this,* based on the Man- 
delstam representation, involves the solution of a 
very complicated set of coupled integral equations 
on a computing machine, and is therefore not very 
transparent. One qualitative feature which has 
been emphasized by Chew and Mandelstam,° how- 
ever, is that the existence of the resonance seems 
to follow from the operation of a “bootstrap mech- 
anism,” in which the strong force between two pi- 
ons in a P state, which is needed to produce the 
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resonance, is provided by the exchange of a pair 
of resonating pions. 

In this way, the basic underlying source of the 
pion-pion interaction (which could, for example, 
be the A¢* interaction or a force produced by the 
exchange of a strongly interacting S-wave pion 
pair) does not seem to play a large quantitative 
role in the final results, but merely provides the 
spark which sets the bootstrap off. One would, 
therefore, expect to be able to obtain the p mes- 
on’s properties as the result of a self-consistent 
calculation with no parameters. 

It is our purpose here to apply the bootstrap 
mechanism in a very simple-minded and trivial 
calculation, which yields quantitative values for 
the p mass and coupling constant. The approxima- 
tion is the following: ap meson, of mass My and 
coupled to the pion with a coupling constant ypz_; 
is exchanged between two pions as shown in Fig. 1. 


FIG, 1. The one p-meson exchange diagram. 
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The force thereby produced in the J=1, T=1 
channel is attractive and depends on mp and yprr. 
If m, and Yp__ are adjusted properly, this force 
may be made to produce a resonance at m, with 
coupling constant ypq7z. Thus we obtain two rela- 
tions between mp and ypr_ from which both may 
be determined. 

In practice, the calculation is most easily car- 
ried out using the determinantal method.” One 
computes the determinantal function for the J=1, 
T=1 channel from the equation® 

2 
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(1) 
where s is the total c.m. energy squared, 
i6 
ImD = -(sinée ),, (2) 


and (sinde*®), is calculated from the graph of Fig. 
1. Then we require 


ReD(m p /™,” Yona? ™ p/m + =0, (3) 
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Equations (3) and (4) together fix mp and yprz. 
The computation is explicit and straightforward, 
and an approximate numerical evaluation of the 
necessary two integrals yields 


2 
m = 950 Mev, ‘oes /4n = 2.8. 








There are no parameters to be adjusted in obtain- 
ing these results, other than the pion mass which 
only provides a dimension. The numbers are in 
fair agreement with the present experimental data,' 
which indicate something like 


2 
m, = 750 Mev, pom /4n = 1. 


The most serious deficiency of a calculation like 


this based on the determinantal method lies in the 
violation of crossing symmetry. One may hope, 
however, in analogy to other more or less success- 
ful calculations,’ that this lack is not of great quan- 
titative importance; in any case, it may be of val- 
ue to illustrate by a simple method what it is that 
the complicated and more complete calculations 
are supposed to do if they can ever be successfully 
carried out. 

In conclusion, it may be remarked that the “boot- 
strap” mechanism may be responsible for the exis- 
tence of many, if not all, other particles as well; 
for example, the strange vector meson called M 
by Gell-Mann,° with isotopic spin 3, may be a P- 
wave resonance in K7 scattering. There may, in 
addition, be an S-wave resonance or pseudoreso- 
nance, corresponding to a K’ meson like the K 
but with opposite parity. In a complete calculation 
of the bootstrap in the 77 system, all of these par- 
ticles, or resonances, should of course be in- 
cluded. For a more complete discussion of these 
points, see reference 5. 
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numerical calculations. 
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